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Preface 


First-degree or, to use the generally accepted term, linear inequal- 
ities are inequalities of the form 


ax +by+c>0 


(for simplicity we have written an inequality in two unknowns x 
and y). The theory of systems of linear inequalities is a small but 
most fascinating branch of mathematics. Interest in it is to a consider- 
able extent due to the beauty of geometrical content, for in geoniet- 
rical terms giving a system of linear inequalities in two or three un- 
knowns means giving a convex polygonal region in the plane or a 
convex polyhedral solid in space, respectively. For example, the 
study of convex polyhedra, a part of geometry as old as the hills, 
turns thereby into one of the chapters of the theory of systems of 
linear inequalities. This theory has also some branches which are 
near the algebraist's heart; for example, they include a remarkable 
analogy between the properties of linear inequalities and those of 
systems of linear equations (everything connected with linear equations 
has been studied for a long time and in much detail). 

Until recently one might think that linear inequalities would for- 
ever remain an object of purely mathematical work. The situation 
has changed radically since the mid 40s of this century when there 
arose a new area of applied mathematics — linear. programming — 
with important applications in the economy and engineering. Linear 
programming is in the end nothing but a part (though a very impor- 
tant one) of the theory of systems of linear inequalities. 

It is exactly the aim of this small book to acquaint the reader 
with the various aspects of the theory of systems of linear inequali- 
ties, viz. with the geometrical aspect of the matter and some of the 
methods for solving systems connected with that aspect, with certain 
purely algebraic properties of the systems. and with questions of 
linear programming. Reading the book will not require any know- 
ledge beyond the school course in mathematics. 

A few words are in order about the history of the questions to 
be elucidated in this book. 

Although by its subject-matter the theory of linear inequalities 
should, one would think, belong to the most basic and elementary 
parts of mathematics, until recently it was studied relatively little. 
From the last years of the last century works began occasionally to 
appear which elucidated some properties of systems of linear inequal- 
ities. In this connection one can mention the names of such mathe- 
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ici i wski (one of the greatest geometers of the 
E A ae of this century especially well 
tien for his works on convex sets and as the creator of “Minkow- 
skian geometry”), G.F. Voronoi (one of the fathers of the : ra 
tersburg school of number theory ), A. Haar (a Hungarian mathe- 
matician who won recognition for his works on “group integration 
H. Weyl (one of the most outstanding mathematicians of the first ha 
H, this century; one can read about his life and work in the pamph- 
let “Herman Weyl” by I. M. Yaglom, Moscow, “Znanie”, 1967). 
Some of the results obtained by them are to some Extent or other ref- 
lected in the present book (though without mentioning the authors 
name not until the 1940s or 1950s, when the rapid growth of applied 
disciplines (linear, convex and other modifications of “mathematical 
programming”, the so-called “theory of games etc.) made an advanced 
and systematic study of linear inequalities a necessity, that a really 
intensive development of the theory of systems of linear inequali- 
ties began. At present a complete list of books and papers on inequal- 
ities would probably contain hundreds of titles, 


1. Some Facts from Analytic Geometry 


T. Operations on points. Consider a plane with a rectangular coordi- 
nate system. The fact that a point M has coordinates x and y in 
this system is written down as follows: 


M — (x, y) or simply M(x, y) 


The presence of a coordinate s 
operations on the points of th 


addition of points and the 
number. 


The addition of points is defined in the following way: if M, = 
(X1, yı) and M3 — (x5, y2), then 


ystem allows one to perform some 
€ plane, namely the operation of 
operation of multiplication of a point by a 


M; + M;-(xi xs, y, +2) 
Thus the addition of points is reduced to the 
similar coordinates. 
The visualization of this operation 
point M; + M; is the fourth vertex of 


ted on the segments OM, and OM, as 
coordinates). M,, 
parallelogram. 


8 


addition of their 


is very simple (Fig. 1); the 
the parallelogram construc- 
its sides (O is the origin of 
O, M; are the three remaining vertices of the 


The same can be said in another way: the point M, +M, is 
obtained by translating the point M; in the direction of the segment 
OM, over a distance equal to the length of the segment. 

The multiplication of the point M(x,y) by an arbitrary number 
k is carried out according to the following rule: 

kM — (kx, ky) 
The visualization of this operation is still simpler than that of the 
addition; for k > 0 the point M' — kM lies on the ray OM, with 


D +H, 
"% 


Ü E 
* Fig. ! 
OM' =k x OM; for k <0 the point M' lies on the extension of 
the ray OM beyond the point O, with OM’ = |k| x OM (Fig. 2). 
The derivation of the above visualization of both operations 
will provide a good exercise for the reader". 


"TH (k<0) 
Fig. 2 


The operations we have introduced are very convenient to use 
in interpreting geometric facts in terms of algebra. We cite some 
examples to show this. 


* Unless the reader is familiar with the fundamentals of vector theory. 
In vector terms our operations are known to mean the following: 
the point M, + M; is the end of the vector OM, + OM; and the point kM 
is the end of the vector k x OM (on condition that the point O is the 
beginning of this vector). 
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(1) The segment M,M; consists of all points of the form 
sıM 1 53M5 


where s,, s; are any two nonnegative numbers the sum of which 
Sy, 5 y twi 

equals 1. T : i 

ere a purely geometric fact, the belonging of a point to the 
segment M,Mz, is written in the form of the algebraic relation 
M — 5M, + s2M; with the above constraints on s}, 83. 


K M 
d e 
4 
2, M H 
4 $ 
Fig. 3 Fig. 4 


To prove the above, Consider an arbitrary point M on the 
segment M,M>. Drawing through M straight lines parallel to OM; 
and OM, we obtain the point N, on the segment OM, and the 
point N, on the segment OM, (Fig. 3). Let 


M3M MiM 
“> MyM, 2 DIM, M, 


the numbers s; and 52 being nonnegative and their sum equalling 1. 
From the similarity of the corresponding triangles we find 

ON, M3M ON, MM 

OM, ~ MyM, ~ OM, MM, ° 


hence 
that when the point M 
rection from M, toward 
he values from 0 to 1. 


which yields N, = s,M,, Nz=52M). But M= N,+N,, 
M=5,M, +5.M). We, finally, remark 
runs along the segment .M,M), in the di 
M2, the number s, runs through all t 
Thus proposition (1) is proved. 

(2) Any point M of the straight line M,M, can be represented as 


(M; +(1 = ÒM, 


where t is a number. 
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| 


Jj 


In fact, if the point M lies on the segment M,M;, then our 
statement follows from that proved above. Let M lie outside of 
the segment M,M;. Then either the point M, lies on the segment 
MM; (as in Fig. 4) or M; lies on the segment MM,. Suppose, 
for example, that the former is the case. Then, from what has 
been proved, 


M,-sM-c(l—sMa (0<s<1) 


Hence 


M=1M,—+=5M,=1M, ( - 0M; 


s 
where t = I/s. Let the case where M; lies on the segment MM, 
be considered by the reader. 

(3) When a parameter s increases from 0 to œ, the point sB 
runs along the ray OB* and the point A + sB is the ray emerging 
from A in the direction of OB. When s decreases from 0 to — co, 
the points sB and A + sB run along the rays that are supplementary 
to those indicated above. To establish this, it is sufficient to look 
at Figs. 5 and 6. 

It follows from proposition (3) that, as s changes from 
— c» to +, the point 4+sB runs along the straight line 
passing through 4 and parallel to OB. 

The operations of addition and multiplication by a number can, 
of course, be performed on points in space as well. In that case, 


* The point B is supposed to be different from the origin of coor- 
dinates O. - 
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by definition, 


My, + Mz — (X1 + X2, yi + y2, zi + Zo) 


| 
kM = (kx, ky, kz) ` | 

All the propositions proved above will obviously be true for | 

a ell. " r X d 
PNG P Aide this section by adopting a convention which will 
later help us formulate many facts more clearly and laconically. 
Namely, if # and £ are some two sets of points (in the plane’ 
or in space), then we shall agree to understand by their "sum 
JC +L a set of all points of the form K+L where K is an 
arbitrary point in # and Lan arbitrary point in Y, 


j er pre 
/ / j 
/ / / 
n / i g 
y Ge 22 137 er 


Fig. 7 Fig. 8 
Special notation h 


J as been employed in mathematics for a long 
lime to denote the 


belonging of a point to a given set; namely, 
€ that a point M belongs to a set .// one 
Symbol e standing for the word “belongs”). 
Set of all points of the form K--L where 


writes Me.// (the 
So J^ is a 
Ke and Le &. 

From the visualization of the addition of points a simple rule 


for the addition of the point sets 4% and Z 
OK over a distance equal to the length 
the latter that will be H+ PL. 


We shall cite some examples. 
1. Let a set Jf^ consist of a 


I2 


2. X and Z are segments (in the plane or in space) not 
parallel to each other (Fig. 9). Then the set # -- is a paral- 
lelogram with sides equal and parallel to # and ¥ (respectively). 
What will result if the segments # and & are parallel? 

3. # is a plane and ¥ is a segment not parallel to it. The 
set 4 +L is a part of space lying between two planes parallel 
to X (Fig. 10). 


i 


tar 


Fig. 9 Fig. 10 


4. Jf and 4^ are circles of radii r; and rz with centres P; 
and P, (respectively) lying in the same plane n. Then # + ^ is 
a circle of radius r; +r2 with the centre at the point P, + P; 
lying in a plane parallel to x (Fig. 11). 

2°. The visualization of equations and inequalities of the first 
degree in two or three unknowns. Consider a first-degree equation 
in two unknowns x and y: 

ax + by+c=0 (1) 
Interpreting x and y as coordinates of a point in the plane, it is 
natural to ask the question: What set is formed in the plane by 
the points whose coordinates satisfy equation (1), or in short what 
set of points is given by equation (1)? 

We shall give the answer though the reader may already know 
it: the set of points given by equation (1) is a straight line in the 
plane. Indeed, if b #0, then equation (1) is reduced to the form 

y=kx+p 
and this equation is known to give a straight line. If, however, 
b=0, then the equation is reduced to the form 
EN 


and gives a straight line parallel to the axis of ordinates. 
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A similar question arises concerning the inequality 
2 
ax +by+ce>0 (2) 
What set of points in the plane is given by inequality (2)? 


Here again the answer is very simple. If b z 0, then the inequali- i 
ty is reduced to one of the following forms 


y2kxtp or y<kx+p 


It is easy to see that the first of these inequalities is satisfied by 


all points lying "above" or on the straight line y= kx + p and the 
second by all points lying “below” or on the line (Fig. 12). If 
however, b —0, then the inequality 


is reduced to one of the 
following forms 


xzh or x<h 


the first of them being satisfied by all points lying to the "right" 
of or on the straight line x — h and the second by all points to 
the “left” of or on the line (Fig. 13). 

Thus equation (1) gives a straight line in the 
and inequality (2) gives one of the 
this line divides the whole 
belong to either of these t 


coordinate plane 
two half-planes into which 
plane (the line itself is considered to 


wo half-planes). 
We now want to solve simil 


i ar questions with regard to the 
equation 

ax + by +ez+d=0 (3) 
and the inequality 

ax+by+ez+d>0 (4) 


of course, here x, J: 2 are interpreted as Coordinates of a point 
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e 


in space. It is not difficult to foresee that the following result will 
be obtained. 

Theorem. Equation (3) gives a plane in space and inequality (4) 
gives one of the two half-spaces into which this plane divides the 
whole space (the plane itself is considered to belong to one of 
these two half-spaces). 

Proof. Of the three numbers a, b, c at least one is different 


^ 
"€. 
io 


Fig. 12 
from zero; let c #0, for example. Then equation (3) is reduced 
to the form 


z=kx+ly+p (5) 


Denote by ¥ the set of all points M (x, y, z) which satisfy (5). 
Our aim is to show that ^ is a plane. 

Find what points in Z belong to the yOz coordinate plane. 
To do this, set x — 0 in (5) to obtain 


z=ly+p (6) 


Thus the intersection of ^ with the yOz plane is the straight line u 
given in the plane by equation (6) (Fig. 14). 4 

Similarly, we shall find that the intersection of Z with the xOz 
plane is the straight line v given in the plane by the equation 


z=kx+p (7) 


Both lines u and v pass through the point P(0. 0, p). 

Denote by x the plane containing the lines u and v. Show 
that x belongs to the set Z. f : 

In order to do this it is sufficient to establish the following 
fact, viz. that a straight line passing through any point Ae v and 
parallel to u belongs to 2. . 

First find a point B such that OB|u. The equation z — ly +p 
gives the straight line u in the yOz plane; hence the equation z = ly 
gives a straight line parallel to u and passing through the origin 
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(it is shown as dotted line in Fig. 14). We can take as B the 
point with the coordinates y — 1, z=! which lies on this line. 
An arbitrary point Ae v has the coordinates x, 0, kx + p. The 
point B we have chosen has the coordinates 0, 1, /. The straight 
line passing through A and parallel to u consists of the points 


A+sB=(x, 0, kx + p t s(0, 1, D) — 


= (x, s, kx + p + sl) 
where s is an arbitrary number (see proposition (3) of section | ) 


Fig. 13 Fig. 14 


It is easy to check that the coordinates of a point A+ sB 
Satisfy equation (5), i.e. that A+sBeY. This proves that the 
plane x belongs wholly to the set 2. 


It remains to make the last Step, to show that 4 coincides 
with T or, in other words, that the set 7 does not contain any 
points outside zm. 

To do this, consider three points: a point M (Xo. yo. zo) 
lying in the plane m, a point M’ (xo, Yo, Zo +£) lying "above" 
the plane x (e 0), and a point M" (xo, yo, Zo — £) lying *below" x 
(Fig. 15). Since Men, we have zo = kx + lvo + p and hence 


Zo +E> kxo + lyo + p 

Zo — € < kxo + lyo + p 
This shows that the coordinates of the point M* satisfy the 
strict. inequality 


zZ>kx+ly+p 

and the coordinates of the point M" satisfy the strict inequality 
Z<kx+ly+p 
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Thereby M’ and M" do not belong to ^ This proves that Z 
coincides with the plane m. In addition, it follows from our 
arguments that the set of all points satisfying the inequality 


ax by+cz+d20 


is one of the two half-planes (either the "upper" or the "lower" one) 
into which the plane x divides the whole space. 


oM Torat) 


Fig. 15 


2. Visualization of Systems of Linear Inequalities 
in Two or Three Unknowns 


Let the following system of inequalities in two unknowns 


x and y be given 


(0) 


ax + bay Cm SO 


In the xOy coordinate plane the first inequality, determines 
a half-plane Hi. the second a half-plane I, etc. If a pair of 
numbers x, y satisfies all the inequalities (1), then the correspond- 
ing point M(x. y) belongs to all half-planes Ty, Il5 25 b. 
simultaneously. In other words, the point M belongs to the 
intersection (common part) of the indicated half-planes. It is easy 
to see that the intersection. of a finite number of half-planes 
is a polygonal region Jf. Figure 16 shows one of the possible 
Tegions YW. The area of the region is shaded along the boundary. 
The inward direction of the strokes serves to indicate on which 
Side of the given straight line the corresponding half-plane lies; 


the same is also indicated by the arrows. 
17 


The region J^ is called the feasible region of the system (1). 


Note from the outset that a feasible region is not always bounded: 
as a result of intersection of several half-planes an unbounded ` 
region may arise, as that of Fig. 17, for example. Having in mind 
the fact that the boundary of'a region X consists of line 
segments (or whole straight lines) we say that % is a polygonal 
feasible region of the system (1) (we remark here that when 


, j 

$i 

d 

vu 

z i } 

| 

Fig. 17 Fig. 18 a | 

à region .# is bounded it is simply called a polygon*) Of course — ' 


à case is possible where there is not a single point simultaneously 
belonging to all half-planes under consideration, i.e. where the 
region % is "empty"; this means that the system (1) is incom- 
patible. One such case is presented in Fig. 18. 

Every feasible region J^ is convex. It will be recalled that 
according to the general definition a ser of points (in the plane | 
or in space) is said to be convex if together with any two 
points A and B it contains the whole segment AB. Figure 19 


illustrates the difference between a convex and a nonconvex set. 
The convexity of a feasible 


region J^ ensues from the very way 

in which it is formed; for it is formed by intersection of several 
half-planes, and each half-plane is a convex set. 

Lest there be any doubt, however, with regard to the convexity 

of J^, we shall prove the following lemma. 


* Here to avoid any misunderstanding we must make a 
tion. The word “polygon” is understood in the school geometry course 
to designate a closed line consisting of line segments, whereas in the litera- 
ture on linear inequalities this term does not designate the line itself but 
all the points of a plane which are Spanned by it (i. e. lie inside or on 


the line itself). It is in the sense accepted in the literature that the term 
"polygon" will be understood in what follows: 
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Lemma. The intersection of any number of convex sets is 
a convex set. 

Proof. Let ^, and -#, be two convex sets and let # be 
their intersection. Consider any two points 4 and B belonging 
io W (Fig. 20). Since 4e2£,, Be, and the set J£, is convex, 
the segment AB belongs to #,;. Similarly, the segment AB 
belongs to .# . Thus the segment AB simultaneously belongs to 
both sets 4%, and -#, and hence to their intersection A. 


Fig. 20 


This proves that J£ is a convex set. Similar arguments show 
that the intersection of any number (not necessarily two) of convex 
Sets is a convex set. 

Thus the locus of the points whose 
inequalities (1), or equivalently the feasible 
is a convex polygonal region X. It is a resu 
all half-planes corresponding to the inequalities of the given system. 

Let us turn to the case involving three unknowns. Now we 


coordinates. satisfy all the 
region of the system (1). 
lt of intersection of 
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are given the system 
ax+by+coz+d,>0 


ax + boy + ez +d, 20 O 


aX + bY +¢,2 +d, >0 


m 


As we know from Section 1, each of the above inequalities 
gives a half-plane. The region determined by the given system 


A E 
4 ETTER] 
Tae nl Fig. 22 Fig. 23 


will therefore be the intersection (common part) of m half-planes, 
and the intersection of a finite number of half-planes is a convex 
polyhedral region 2. Figure 21 exemplifies such a region for m — 4. 
In this example the region .% is an ordinary tetrahedron (more 
strictly, # consists of all points lying inside and on the boundary 
of the tetrahedron); and in general it is not difficult to see that 
any convex polyhedron’ can be obtained as a result of the 
intersection of a finite number of half-planes.* Of COUrse, a case 
is also possible where the region J^ is unbounded (where it 
extends into infinity); an example of such a region is represented 
in Fig. 22. Finally, it may happen that there are no points 
at all which satisfy all the inequalities under consideration (the 
System (2) is incompatible): then the region J^ is empty. Such à 
case is represented in Fig. 23. 


| of all points of space spanned by 
ace itself. the set of course including 
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Particular attention should be given to the case where the 
system (2) contains among others the following two inequalities: 


ax + by+cz+d>0 

—ax—by—cez—dz0 

Which can be replaced by the sole equation 
ax+by+cz+d=0 


The latter gives a plane m in space. The remaining inequalities (2) 


will separate in the plane x a convex polygonal region which 
will be the feasible region of the system (2. It is seen that 
a particular case of a convex polygonal region in space may be 
represented by a convex polygonal region in the plane. In Fig. 24 
the region of is a triangle formed by the intersection of five 
half-planes, two of them being bounded by the "horizontal" 
plane m and the remaining three forming a "vertical" trihedral 
prism, 

By analogy with the case involving two unknowns, we call 
the region # the feasible region of the system (2). We shall 
emphasize once again the fact that a region # being the 
intersection of a number of half-planes is necessarily convex. 

Thus the system (2) gives a convex polyhedral region X in space. 
This region results from intersection of all half-planes corresponding 
to the inequalities of the given system. 

If the region 2 is bounded, it is simply called the feasible 
polyhedron of the system (1). 


3. The Convex Hull 
of a System of Points 


Imagine a plane in the shape of an infinite sheet of plywood to 
have pegs driven into it at points A1, Az, ..., Ap Having made 
a rubber loop, stretch it well and let it span all the pegs (see the 
dotted line in Fig. 25). Then allow it to tighten, to the extent 
permitted by the pegs of course. The set of the points spanned by 
the loop after the tightening is represented in Fig. 25 by the shaded 


Fig. 26 


area. It appears to be a convex polygon. The latter is called the 

convex hull of the system of the points Ay, Az, n Ap. 

If the points A,, Ap, ..., A, are located in space rather than in the 
plane, then a similar experiment is rather hard to put into practice. 
Let us give rein to our imagination, however, and suppose that we 
have managed to confine the points Aj A, ., A, ina bag made 
of tight rubber film. Left to its own devices, the bag will tighten 

, until prevented from doing so by some of the points. Finally, a time 
will arrive when any further tightening is no longer possible (Fig. 20) 
It is fairly clear that by that time the bag will have taken shape 
ofa convex polyhedron with vertices at some of the points A4, Az, -s Ap: 
The region of space spanned by this polyhedron is again called the 
convex hull of the system of the points Ai, A5, Age 

Very visual as it is, the above definition of the convex hull is no! 
quite faultless from the standpoint of “mathematical strictness”. We - 
now define that notion strictly. 


Let 4i, A>, ., Ap be an arbitrary set of points (in the plane 
or in space). Consider all sorts of points of the form 


8141 F $545 + ..4 SpA a) 


N 
N 


ae e 


that Ae., So any of the above segments belongs wholly to M, 


where sy, S2, ..., Sp are any nonnegative numbers whose sum is one: 


Sia 52, «3 Spz^ 0! and Si 4232-5 at Spe (2) 
Definition. The set of points of the form (1) with condition (2) is 
called the convex hull of the system of the points Ai, Az, -= Ap 


and denoted by 
(Av, Az, s Ap? 


To make sure that this definition does not diverge from the former, 
first consider the cases p — 2 and p — 3. If p — 2, then we are given 
two points A; and 45. The set (4, 42), as indicated by proposition 
(1) of section 1, is a segment 4,45. 

If p — 3, then we are given three points 4;. A2 and 43. We show 
that the set (Aj. 42, A3) consists of all the points lying inside and 


on the sides of the triangle 444243. 


Moreover, we prove the following lemma. 

Lemma. The set (Aj, n Ap 1, Ap? consists of all segments joining 
the point A, with the points of the set Ai, = Ap-1>: 

Proof. For notational convenience, denote the set <A1, «s Ap-12 
by -%-1 and the set (A, n Api. Ap? by Mp 

Consider any point Ae./%,. It is of the form 


A= SA, +. + Sp-1AÀp-1 + Spâp 


where 
Sin xe Spa Ur Si en S UE 1 


If 5, = 0, then A €.7,.,: thus the set ./,. , is a part of p If s, — 1. 


then A = A,; thus the point A, belongs to ./£. So M, contains ./f, .. , 
and the point 4,. We now show that any segment A’A, where 


A'e.f/,_, belongs wholly to . 4. 
If A is a point of such a segment, then 


A=tdA'+sA, (t s>0,t+s=1N 
On the other hand, by the definition of the point A’ we have 
Al = tyAy + + tp-14p-1 
Deen a ia we + tpa = 
hence 
A= ttjAy ES + tip Aa SA, 


we have (1). (2). This proves 
(4 


Setting tty = St. es Hp = Sp S — Sp 


29) 


where 


Now it remains for us to check that the set .//, does not contain 
anything but these segments, i.e. that any point of ./£, belongs to 
one of the segments under consideration. . 

Let Ac.//,. Then we have (1), (2). It can be considered that s, # l, 
otherwise A = A, and there is nothing to be proved. But if s, # 1. 


then s; +..+5,-,;=1—s,>0, therefore we can write down 
2 4, + 
A — (s FASES a) aay jobs 
opal A +s,A 
pr P 
eee w= pur 


The expression in square brackets determines some point A 
belonging to .//, _ , for the coefficients of Aj, .., 4,,.. , are nonnegative 
in this expression and their sum is one. So 


A — (81 5-1) A + SpA, 


Since the coefficients of A’ and A, are also nonnegative and their 
sum is one, the point A lies on the segment A’A,. This completes 
the proof of the lemma. 

Now it is not difficult to see that the visual definition of the 
convex hull given at the beginning of this section and the strict 
definition which follows it are equivalent. Indeed, whichever of the 
two definitions of the convex hull may be assumed as the basis, 
in either case going over from the convex hull of the system A), ~ 
A,-1 to that of the system Aj, ..., A,~1, A, follows one and the 


same rule, namely the point A, must be joined by segments to 
all the points of the convex hull for Aj 2. Ap-, (this rule is 


immediately apparent when the convex hull is visually defined and 
in the strict definition it makes the content of the lemma). If we 
now take into account the fact that according to both definitions 
we have for p — 2 one and the same set, the segment 4,45, the 
equivalence of both definitions becomes apparent. 

The term *convex hull" has not yet been quite justified by us, 


however, for we have not yet shown that the set (A,, Ag, ws Ap? 
is always convex. We shall do it now. 


Let A and B be two arbitrary points of this set: 
A-— 541-4 $5345 +... + SpA, 
B=t,A,;+t,4,4+...4 UA, 


Si Sp Lay ony ty > 0 


to 


Pm 


> 


Any point C of the segment AB is of the form 
C=sA+tB 
where 
s,t>0, str=1 (4) 
which gives r 
C — s(sqAi +. + SpAp) + tliA +. + (AP) — 
= (ss, + ttj) A, +... + (ssp + tty) Ap 


The numbers preceding A, .... Ap as coefficients are nonnegative 
and their sum amounts to one (which follows from (3), (4)). This 
means that the point C belongs to the set (Ai. Az, ~ Ap, i.e: 
this set is convex. 

At the same time it is easy to see that (Ay, A2, ~» Aj? is the 
least of all convex sets that contain the starting points Aj, A, om Ap 
viz. that it is contained in any of these sets. This statement follows 
directly from the lemma proved above and from the definition of 
a convex set. 

This explains the name “convex hull” and at the same time 
provides yet another explanation of the fact why the set (Aj, A2, ~. 
Ap? can be obtained using the method described at the beginning of 
the section. For the set spanned by a rubber loop (or film) after 
the latter has tightened to its limit around the system of the points 
Ais Hay xy Ap S exactly the least convex set that contains the 
indicated points. 


4. A Convex Polyhedral Cone 


Let’s begin with a definition. 

A convex polyhedral cone is the intersection of a finite number 
of half-spaces whose boundary surfaces pass through a common point: 
the latter being called the vertex of the cone. | 

We shall first of all point out how the notion of a convex 
polyhedral cone is related to systems of linear inequalities. We shall 
confine ourselves to a particular case, namely a case where the 
vertex of a cone is the origin of coordinates. This means that all 
boundary planes pass through the origin, and the equation of a 
plane passing through the origin is of the form 


ax + by +cz=0 


(the absolute term in the equation must be equal to zero, otherwise 
(0, 0, 0) will not satisfy the equation). Thus a convex po- 
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lyhedral cone with the vertex at the origin is the feasible region 
of a system of homogeneous inequalities: 


AmX + by + cuz Z0 


The reverse is true of course; the feasible region of a homogencous 


Fig. 27 


system of inequalities is always a convex polyhedral cone with the 
vertex at the origin. 

An example of a convex polyhedral cone may be provided by 
à convex region in space whose boundary is a polyhedral angle 
with the vertex S, a kind of infinite convex pyramid without a base 
and extending unboundedly from the vertex (Fig. 27 represents one 
of such pyramids which has four faces). Other, less interesting cases 
are possible; for example: 

1. A half-space (Fig. 28. a). In such a “cone” the role of the 
vertex may be played by any point S e x, where z is the boundary 
plane of a given half-space. 

2. The intersection of two half-spaces whose boundary planes 
intersect along a straight line / (Fig. 28, b). The role of the verte 
may be played by any point Sel. 

3. A plane. It is clear that any plane m in space can be 
considered as the intersection of two half-spaces lying on different 
sides of z (Fig. 28, c). In this case the role of the vertex may 
be played by any point S e q. 

4. A half-plane (Fig. 28, d). The vertex S is any point of the 
boundary line. 
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_ 3. A straight line. Every line / in space can be obtained by 
intersection of three half-spaces whose boundary planes pass 
through | (Fig. 28, e). The vertex S is any point of the line l. 


@ p e E 
NET, 
(e) Y. 
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J 
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6. An angle (less than 180°) in an arbitrary plane x (Fig. 28, /). 
One can obtain an angle by intersecting the plane 7 with two 
half-spaces (precisely how?). 

7. A ray (Fig. 28, g). A ray can be considered as the intersection 
of a straight line and a half-space. The vertex S is the beginning 
of the ray. 

8. A point. This “cone” can be obtained by taking the common 
part of a ray and the corresponding half-space (Fig. 28, ^). 

Of course, the enumerated examples 1-8 diverge (to a greater 
Or lesser extent) from the usage of the word “cone”, but we are 
compelled to reconcile ourselves to this if we are to preserve the 
general definition of a convex polyhedral cone given at the begin- 
Ning of this section. 

We now try to show in a few words that the sets pointed out 
above exhaust all polyhedral convex cones in space. 
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Let p denote the number of half-spaces whose intersection is the 
cone # under consideration. If p — 1, then our statement is valid. 
for then .% is a half-space. A simple argument which is left to the 
reader shows that, if our statement is valid for a cone which 
results from intersection of p half-spaces, it is also valid for a cone 
formed by intersection of p+ 1 half-spaces. It follows according 
to the principle of complete mathematical induction that our state- 
ment is valid for any p. 

Convex polyhedral cones possess many interesting properties. It 
is beyond the scope of the present book to go into these 
subjects, so we shall confine ourselves to a few simplest propositions. 

We introduce one more definition or notation, if you please. 

Let Bi, Bj .., B, be an arbitrary set of a finite number 
of points (in space). The symbol (B,, B5, . B,) will denote a set 
of points of the form 


tB, + t2B,+..4+1,B 


gq 
where ty, t2, ~s t} are arbitrary nonnegative numbers. 

What is the geometric meaning of the set (B, Bz, .. Bj? 
It is clear from the definition that it is the sum of the sets 
(Bi) (B5). -= (B); we must first see, therefore, what is the 
geometric meaning of the set (B), ie. of the set of points 
of the form 1B, where t is any nonnegative number and B a fixed 
point. But the answer to the last question is obvious: if B 
is the origin, then the set (B) also coincides with the origin; 
otherwise (B) is a ray emerging from the origin and passing through 
the point B. Now we remark that the sum of any set and the origin 
is again the same set; hence it is clear that when studying the sets (Bi: 
Ba, .., B) we shall lose nothing if we consider all the points 
By, B>, ..., B, to be different from the origin. Then the set (B,, Bz 
- B,) will be the sum of the rays (By), (B3), ... (B,). 

The last remark makes the following lemma almost obvious. 

3 Lemma. The set (B,, Bie 
joining each point of the set (B,, -. B, 4) with each point of 
the ray (B,). j E 

The strict proof of the lemma is carried out according to the 
same plan as the proof of the similar lemma of Section 3: 
the reader is advised to carry it out independently. 

It is easily deduced from the lemma that (B,, B5) is an angle. 
a straight line or'a ray (Fig. 29, a, b, c). It is then readily 
established that (Bi, B5, B3) is one of the following sets: 8n 
infinite trihedral pyramid, a plane, a half-plane, an angle, a stra- 
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B,) is a union of segments 


ight line or a ray. Now it becomes clear that there must exist 
a close relation between the sets (Bi) (B2), .. (B,) and convex | 
polyhedral cones. Such a relation does in fact exist. For greater 
intelligibility we shall formulate the corresponding propositions as 
two theorems. 

Theorem 1. The set (Bi, B5, .... B,) either coincides with the whole 
of space or else is a convex polyhedral cone with vertex at 
the origin. 


2 

4 ^ 

E. e 5, 

o Yy 4 0 4 0 
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That the set (Bı, B, ... B,) can indeed coincide with the whole 
of space is shown by the following example. Consider four points 
B,, B, Bs, B, lying in such a way that the rays (B4). (B2), (B3), (Ba) 
form pairwise obtuse angles (Fig. 30). Each of the sets (Bi, B2, Bs), 
(Bı, Bj, B4), (Bi, B3» Ba), (Bi, Bi, Bs) is an infinite trihedral 
pyramid with vertex at the origin. The set (B;, B2, Bs, B4) apparently 
contains all of these pyramids; and the union of the pyramids 
coincides with the whole of space! 
Theorem 2. Any-convex polyhedral cone with vertex at the origin 
is a set of the form (Bi, Bz, .., Bj). 
The proof of Theorem 1 will be carried out on general lines. 
i We shall take advantage of the method of complete mathematical 
induction. The statement of the theorem for q = 1 is obvious. Now 
; we suppose that the theorem is valid for sets of the form (By, .... B,) 
and, guided by this fact, prove it to be valid for the sets 
) (Be Bs Bira): $ 
According to the induction hypothesis, (By. -> B) is either 
the whole of space or a convex polyhedral cone in it. As to the 
first case there is as a matter of fact nothing to be proved in it, 
for then (Bi, .., B, B,,,) is also the whole of space. Let the 
second case occur and (B, ., B, be a convex polyhedral 
cone Jf. According to the lemma, the set (Bi, ~» B, By) 
is a union of segments joining each point of the set KA with 
each point of the ray (B,,,); and, as shown earlier, any convex 


29 


JG" wu ^ —— 


polyhedral cone J^ is either an infinite convex pyramid or one of 

- the sets 1 to 8. Having considered the above union of segments 
for each of these cases, it is not difficult to make sure (check 
it for yourself) that either it coincides with the whole of space 
or is again a convex polyhedral cone. Thus the theorem is true 
for sets of the kind (B,) and for (Bj, ..., B,» B,,,) as well, as soon 
as we suppose it to be valid for (B. ... B). Hence it follows that 
the theorem is true for any q. 


Fig. 30 Fig. 31 

Proof of Theorem 2. Let J^ be a convex polyhedral cone with 
vertex at the origin O. As already stated, .# is either an infinite 
convex pyramid or one of the sets | to 8. 

Let .# be a pyramid. We choose a point on each of its edges 
to get a system of points By, Bj, .. B, We state that the 
set (By, Bz, .... B,) is exactly Af. 

To prove this, consider a plane x intersecting all the edges 
of the pyramid J/. We get the points By, Bj, ., B' (Fig. 31). 
Apparently. : 3 

B| = kiBis Bi—kjB B, = k;B, (1) 


where ky, kz, .., k, are nonnegative numbers. 


Now suppose B is a point of the pyramid different from the 
vertex O. The ray OB intersects with the plane z at a point B'. 


It is obvious that B’ belongs to the convex hull of the system 
Bj. Bj, ng B, and hence 


B’ = s,Bí +s2B3 +... + sB; 


WHere.Sys, S2, c9 s, are nonnegative numbers whose sum is One- 
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Now taking into account (1) we have 
B’ = s,k,B, + s2k2B2 +...4+ s,k;B, 


and if in addition it is remembered that B'— KB (k» 0) we 
find that 


B = t,B; + t5B; + —. + tB, 


where 


Gay (Ï=1, 2, 9) 


Thus we have shown that any point B of the pyramid # belongs 
to the set (Bj, Ba, ..., B,). The converse (i.e. that any point of the 
set (Bj, Bz, ., B) belongs to #) is obvious. So J^ coincides 
| with (Bi, Ba, e B). 
The case where .# is one of the exceptional sets 1 to 8 can be 
h proved without much trouble and is left to the reader. 


5. The Feasible Region of a System 
of Linear Inequalities in Twv Unknowns 


Our task now is to give an effective description of all solutions 
of a system of linear inequalities. In the present section we shall 
deal with systems involving two unknowns x and y. In spite of 
the fact that the number of unknowns is not large (only two), 
we shall try to carry out the analysis of these systems from 
general positions, i.e, so that the resülts obtained may be easily 
extended to systems in a larger number of unknowns. 

The solution of any system of linear inequalities is in the long 
run reduced to the solution of a number of systems of linear equations. 
We shall regard the solution of a system of linear equations 
as something simple, as an elementary operation, and shall not be 
confused if, to realize the proposed method, we have to perform 


" this operation many times. * 4 
1°. The necessary lemmas. Given a system. of inequalities 


ax by ei 20 
ü4x + by - c5 20 
| CMM () 
M o ax + b,y + c, 2:0 


It is found advisable to consider side by side with it the 
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corresponding system of homogeneous inequalities 
ax +bıy>0 
ax + boy >0 


(2) 


ax Day SO 


as well as the corresponding system of homogeneous equations 
ax +bıy=0 
axx +bay=0 ) 


dX + by = 0 
We denote the feasible region of the system (1) in the xOy 
coordinate plane by .#, that of the system (2) by Xo and that of the 
system (3) by X. Obviously, Z c Jo, where the symbol c stands for 
the words "is contained in"*. 
Lemma 1. The following inclusion holds 
H+ eK 
i.e. the sum of any solution of a given system of inequalities and any 
solution of the corresponding homogeneous system of inequalities is 
again a solution of the given system. 


Proof. Let A be an arbitrary point of J^ and let B be an 
arbitrary point of Xo. Then the following inequalities are valid 


axa + biya +e, z0 a,Xg + biy > 0 
a2X4 + bry, c5 >0 aX, + bry, >0 
and 


AnX, + b,y, + ©, 


Adding each inequality written on the left to the corresponding 
inequality on the right we have 

y(X4 + Xp) + biya + y) +e, > 0 

aXx, + xy) + bia + yg) +c. > 0 


m 20 anXp buys > 0 


* One should not confuse the symbol c with the symbol e introduced 
earlier. The latter is used when speaking of the belonging of a point to a set. 


If, however, one wants to record the fact that one set is a part of another, 
the symbol c is used. 
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These inequalities imply that the pair of numbers X4 + Xg, V4 + Vp, 
the coordinates of the point A + B, is a solution of the original 
system (1), i.e. that A + Be.#. Thus the lemma is proved. 

Lemma 2. (1) If a ray with the beginning at the point A belongs 
wholly to the set A and P is an arbitrary point of the ray, 
then P — Aes. 

(2) If a straight line belongs wholly to # and A, P are 
two arbitrary points of the line, then P — Ae £P. 


P 
I 
A i 
1 
I I 
I I 
] 4 
I 
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Fig. 32 


Proof of (1) Denote the point P — A by B. The ray under 
consideration consists of points of the form 


A+sB (4) 
where s is an arbitrary nonnegative number (Fig. 32). Any of these 
points is, according to the premises, a solution of the system (1), i.e. 

a,(x4 + 3X5) + bilya sy) + ci 20 
(x4 + sxg) + b(y4 + Syp) + c2 > 0 6 


A(X, + SXB) + by, + SYB) + Cn 2 0 


m 


Consider, for example, the first of these inequalities. It can be 
written down in the following form 


(aixa + biya F c1) + s(ajxy + bia) > O 


Since this inequality holds for any s > 0, the coefficient of s must 
clearly be a nonnegative number: 


dX + biyg> 0 
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Similarly, consideration of the remaining inequalities (5) leads to 
d3Xg + boy, > 0 


üsXg + buys SO 


whence we see that the point B belongs to the set Jf^,. 

The proof of (2) is carried out in a similar way. The straight line 
under consideration consists of points of the form (4) where s is 
an arbitrary number. Therefore inequalities (5) are valid for any 
value of s. Hence it follows that in each of these inequalities the 
total coefficient of s must be equal to zero, i.e. 


Xp + biy, — 0 
d3Xg + bayg = 0 
Ay Xp + b,yg = 0 


Therefore Be ^. Thus the lemma is proved. 

It is easy to see that Lemmas | and 2 are valid for systems 
involving any number of unknowns. 

2*. The case where the system of inequalities (1) is normal. Consider 
again the system of inequalities (1) and the corresponding system of 
homogeneous equations (3) The latter system has an obvious 
solution x — 0, y —0 which is called a zero solution. In order to 
investigate the system (1) it turns out to be important to know if 


the system (3) has any nonzero solutions either. In view of this 
we introduce the following 


Definition. A system of linear inequalities is said to be normal 
if the corresponding system of linear homogeneous equations has only 
a zero solution. 

In other words, a system of inequalities is normal if the set ^. 
the feasible region of the corresponding homogeneous system of 


equations, defined above contains only a single point (the origin 
of coordinates). 


The concept of normal system is meaningful with any number of 
unknowns, of course. 


It is not difficult to show that a compatible system of inequali- 
ties is normal if and only if its feasible region X contains no 
straight lines. 

Indeed, if the system is normal, ie. the set Y contains only 
the origin of coordinates, then the region # does not contain any 
Straight lines, which follows directly from the second statement 
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of Lemma 2. If the system is not normal, then the set & 
contains at least one point B different from the origin of coordina- 
tes. Of course, all points of the form kB, where k is any number, 
also belong to #*. But in this case, whatever the point Pe.” 
(and such a point is sure to exist for the system is compatible 
and the region J£ is therefore not empty), the set of all points of the 
form P 4- kB (where k is any number) belongs, according to Lemma 
1, to X. We know that the set is a straight line. So when 


a system is not normal, the region .# contains à straight line. 
This completes the proof of the above proposition. 

In this section we shall study the feasible region of the system 
(1) under the supposition that the system is compatible (the re- 
gion J^ is not empty) and normal. 

From the fact that the region Jf^ does not contain any straight lines 
it first of all follows that it must necessarily have vertices. The 
term “vertex” is understood by us in the following sense (close 
to the intuitive understanding of the word "vertex"). 

The vertex of a region KX is a point of the region which is 
not an interior point for any segment which lies wholly in A. In 
other words, the vertex is a point Ae.# which has the following 
property: any segment belonging to J£ and passing through A must 
have its beginning or end at this point (Fig. 33, a and b, where 
the point A is one of the vertices: in Fig. 33, b the region 
# is a segment). 

We now explain at greater length why the convex set Jf we are 
interested in has vertices. If % lies on a straight line, then it is 
either a separate point, a segment or a ray, and the existence 
of vertices is evident. If ^ does not lie on a straight line, however, 


* Jf the numbers x, y, =. the coordinates of the point B. satisfy a 
homogeneous system. of equations, then the numbers kx. ky. kz, the 
coordinates of the point KB, satisfy this system too. 
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consider the boundary of this set. It consists of segments and rays 
(# does not contain any whole lines). The end of any of these 
segments and the beginning of any of the rays will evidently be the 
vertices of 2. 

One can find the vertices of the region 2£^ without much trouble. 
Notice first of all that in the xOy coordinate plane the ith inequality 
of the system (1) corresponds to the half-plane whose boundary 
line l; is given by the equation 

ax+by+ce,=0 (i=1, 2, .., m) 


Evidently, the point A of the region J^ is a vertex if and only 
if it belongs to two different boundary lines. 
Let us agree to call regular any subsystem .of two equations of 
the system 
ax bye —0 
ax + by 4-05; — 0 
2 2y (6) 


OX + by Cy = 0 


provided it has a unique solution (x,y). 

From the above description of the vertex now results the following 
method of finding the vertices of the region X. 

In order to find all the vertices one should find the solutions 
of all regular subsystems of the system (6) and pick out those 
which satisfy the original system (1). 3 

Since the number of regular subsystems does not exceed C;, the 
number of combinations of m things 2 at a time, the number of 
the vertices of the region # cannot be greater than that either. 
So the. number of the vertices is finite. 


Remark. It follows from the above that, if the region 4% of the 
solutions of a normal system has no vertex, it is empty — the 
system has no solutions (it is incompatible). 

Example 1. Find all the vertices of the region .# given by the 
system of inequalities 


x+ y+120 

x-—2y—-220 

2x— y—-420 
On solving the subsystems 
x+ yti-0| x+y+1=0 PM 
Pom y Me 2x — y-420 
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(all of them prove regular) we find three points: 


E (0.—1) (1; —2) (2,0) 
of which only the second and the third satisfy all the inequalities 
given. So it is the points 
A,(1, —2) and 4A;(2, 0) 


that are the vertices of the region X. 
Let us return to the system (1). Let 


Ai A», .., A 


be all the vertices of the region J£. The set XAy, 4,5 A,>, the 
convex hull of the system of points A,, 45... A,, also belongs 
to AK (for # is a convex region!) But in that case, according 
to Lemma |, the set 


is Abs ac AD EG 


belongs to ^ too. We shall show that as a matter of fact this 
sum coincides with J£, that we have the following 
Theorem. If the system of inequalities is normal, then 


M = Aj, Ay s AD Mo (7) 


where Ay, A... A, are all the vertices of the region Jf. 

Proof. Let P be an arbitrary point of the region. .# different 
from the vertices of the region. The line A,P intersects the convex 
region -# either along a segment 4,4 (Fig. 34) or along a ray 
with origin at A, (Fig. 35) In the second case P— A,6.%o 
(Lemma 2) hence PeA, Jf. In the first case we reason as 
follows, however. If the point A lies on the bounded edge 4,4, of 
the region J^ (as in Fig. 34), then P belongs to the convex hull 
of the points 4,, A, 4; if, however, the point A lies on an 


P 


unbounded edge with origin at vertex A, (Fig. 36), then, according 
to Lemma 1, we have Ae€4; 4- J£, and thereby Pe(A,, Aj» + 
Ay. Thus in all cases the point P proves to belong to the 
set Au ss Aj» +4 o. Thus the theorem is proved. 
Since we are already familiar with the method of finding vertices, the 
only thing we lack for a complete description of the region X^ is the 
ability to find the region Jy. The latter is the feasible region of the 
homogeneous normal system (2) which we now proceed to describe. 

3°. The homogeneous normal system of inequalities (2). Each of the 
inequalities (2) determines a half-plane whose boundary line passes 
through the origin of coordinates. It is the common part of these 
half-planes that is Jo. 
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Among the boundary lines there are in this case at least two 
different ones (for the system (2) is normal!), Hence A either 
coincides with the origin of coordinates (x — 0. y — 0), or is a ray 
with vertex at the origin of coordinates or an angle smaller 
than 180° with vertex at origin. If-we know two points B, and 
B, lying on different sides of the angle (Fig. 37) then all the 


Fig. 36 Fig. 37 


points of the angle may be written down in the form 
B= tB, + tB: (8) 


pT ti and ta are arbitrary nonnegative numbers: and to find 
BA a and B is not difficult at all, if it is remembered 
at each of them (i) belongs to o, ie. satisfies the system (2). 
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and (ii) lies on the boundary of Xo. i.e. satisfies 
l 1 > Le one of a- 
tions (3). If 4% is a ray, then instead of (8) we ye eeann 


B= tB] (9) 


where Bi is any point of the ray (different from the origin) and z is 
an arbitrary nonnegative number. j4 


Fig. 38 


Example 2. Find the region .# 9 of the solutions of the system 


x 20 
x—2yz0 (10) 
2x — yz0 


as well as the region .% of the solutions of the system in Example 1. 
Solution. The system (10) is normal: for the unique solution of 
the corresponding homogeneous system of equations 


x+ y=0 
x—2y=0 (11) 
2x- y=0 


is (0, 0). 

Choose a point satisfying the first equation of (11) (but different 
from (0, 0)), the point C(— 1, 1), for example. We make sure by a 
simple check that the point C does not satisfy all the inequalities 
(10). therefore, neither the point itself nor any point of the ray OC 
(different from the origin O) belongs to J£. On considering the 
point — C (ie. the point (1, — 1) we find that it belongs to 
K 9. So B, = (1. — 1). The second equation is satisfied by the point 
(2. 1); it is also a solution of the system (10), so that B5 = (2, 1). 
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The region Xo consists of the points (Fig. 38) 
tB, +B. = h(l, — 1) tQ, 1) = (f + 2ta — +h) 


where t, and t, are arbitrary nonnegative numbers. J 
Turning to the system of inequalities of Example 1 we notice 
that it is (10) that is the corresponding homogeneous system. From 


Fig. 39 Fig. 40 
the theorem proved above we;have 
X =A, A2) M, 


where A,(1, —2) and A3(2, 0) are the vertices of the region W: 
So J^ consists of the points (Fig. 39) 


s(1, —2) + (1 — s)(2, 0) + (t4 -2t3, — ti + t3) = 
=(2—s +t, + 2t, —2s — t + (2) 


where s is any number in the interval [0, 1] and 1, t; are any 
nonnegative numbers. 


Example 3. Find the feasible region of the system 


2X— yz0 
—4x *2y20 
x+ yz0 


Proceeding as in Example 2 we find only a single ray: 
B-—t(l, 2)= (t, 2t) (tz 0) 


(Fig. 40). 
Example 4. Find the feasible region of the system 
2x —y 20 
x+y20 
—3x+y20 
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In this case none of the equations 


2x— y= 
x+y=0 
—3x+y=0 


has solutions (except (0, 0)) which would satisfy all the given in- 
equalities. The region % consists of a single point (0, 0), the 
origin of coordinates. 


Fig. 41 


4^. The case where the system of inequalities (1) is not normal. 
This means that the feasible region 7 of the homogeneous system 
of equations (3) contains something else besides the origin of coordi- 
nates. Therefore, all the equations (3) determine the same straight 
line in the plane, and that line is Z. 

According to Lemma 1 the region # contains together with each 
of its points P a straight line P4- (a line passing through P 
and parallel to 7). Consider a straight line .7 not parallel 
to 2. If we know which points of the line .7 belong to 
the region .# — denote the set of these points by J£ ; —then we 
shall be able to find the region .% itself, for then J£ =% +2 
(Fig. 41). 

The equation of the straight line ^ is a,x t biy 0. One of 
the coefficients a, or b, of this equation is different from zero: 
let b, #0, for example. One can then take as a straight line 7 not 
parallel to Z^ the y-axis (its equation is x — 0). In this case the set 
H , —now denote it by J^, —is the part of the y-axis contained 
by #. To find this set, one should put x=0 in the system (1). 


4l 


Then we have the system of inequalities 
byte, 20 
bsy+c.20 (12) 


bay + cp 2 0 


in one unknown y which is easy to solve *. Notice that the set AP, may 
be either an empty set (then .# is also empty) or a point, a segment Or 
a ray (but it cannot be the whole of the y-axis, for otherwise J^ would 
be the whole of the plane, which is impossible). On finding this set we 
shall know the region .¥ itself, for 


HaK, +L (13) 


(if Z is not parallel to the y-axis). 
Example 5. Find the feasible region HA of the system 


^s arbe 


It is easily seen that the system is not normal and that Z^ is the 
straight line 


x+y=0 
(not parallel to the y-axis). Setting x = 0. we have the system 
y-120 
—y+220 
2y+320 
from which it is seen that -#,. the intersection of J£ with the 


y-axis, is a segment with the ends C,(0, I) and C,(0, 2). So He 
is a set of points of the form (Fig. 42) 


(0, y) 4E — x) (x, y — x) 


where x is arbitrary and y is any number in the interval from 
I to 2. 


. * Notice that the system (12) (viewed as a system of inequalities 
in one unknown) is now normal. Indeed, if otherwise, the corresponding 
homogeneous system would have a nonzero solution y*; but then the sys- 
tem (3) would have the solution (0, y*) not belonging to 4. 
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Tn conclusion we shall briefly discuss a theorem which follows from 
the results obtained above. In the two-dimensional case we are consi- 
dering (i.e. where everything takes place in the plane) this theorem is 
not particularly striking and it would be right to regard it as a 
starting point of the extension to the “n-dimensional” case to be 
studied in Section 7. 


Fig. 42 


Theorem. Any (nonempty) convex polygonal region A in the plane 
can be represented as the sum 


4A, v A E DOE (14) 


The augend of the sum is the convex hull of a system of points 4. 
Ax, m» Ap, and the addend is a set of points of the form 
tıBı + Bz +... + GBG where t3, ta, ..., Iq are arbitrary nonnegative 
numbers. 

‘The theorem can be proved in a few words. Consider a system 
of inequalities giving J£. If the system is normal, then equation 
(7) holds: bearing in mind that in this equation Jf is one of the 
sets of the form (Bı, B2), (Bı) or (O) (the origin of coordinates), 
we find that in the case where the system is normal our statement 
is valid. If the system is not normal, then equation (13) holds from 
which it also follows that -4^ is representable in the required form 
(why?) 

Notice that if all points 41. 45... A, coincide with the origin O. 
then the set CA, Az. 4,» also coincides with 0: then only the 
addend. is left of the sum (14). When the points Bi, Bj... B, 
coincide with O. the set (Bj, Bz... Bj) also coincides with O and 
only the augend is left of the sum (14). 
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The converse theorem also holds, though with some reservation. 
Theorem. Any set of the form 


CAT As, a5 Ap) + (Bo Ba, x Ba) 


in the plane is either the whole plane or some convex polygonal region 
in it. À 
The proof is fairly obvious. The addend, i.e. the region Jo = 
(Bi. Bs... Bj). is either the whole plane, or a half-plane, or an angle 
(smaller than 180°), or a ray, or a point (the origin of coordinates). 


whereas the summand ^, = A, Aj... Aj? represents a convex 
polygon. The set J^, +.% can be obtained by translating #o to 
the segments OK, (where K, is any point of .%,) and taking the 
union of the sets obtained (Fig. 43). It is easily seen that this produces 
either the whole plane (this is the case if #5 is the whole plane) 
or some convex polygonal region in it. 


6. The Feasible Region of a System 
in Three Unknowns 


After the detailed analysis given in the preceding section we are 
now in a position to minimize the required theory when considering 
Systems in three unknowns. 

Along with the original system 


ax+by+ez+d,>0 
NOLO DENEN Porra. (1) 
dX buys euz4- d, > 0 


m 
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as in Section 5, we again consider two more systems: 
4ax+by+ez>0 


Ee aona Liebe s CNN I (2) 
Qux + buy + Cuz 


V 
o 


and 
ax+by+cez=0 


per NN (3) 
AX + bay + ez = 0 


We shall again denote the feasible region of the system (1) by %, 
that of the system (2) by # and that of the system (3) by 
2. To use the previously introduced terminology, J^ is some con- 
vex polyhedral region in space and J^; is a convex polyhedral cone. 
As already noted, Lemmas 1 and 2 of Section 5 remain valid here 
too. 

1°. The case where the system of inequalities (1) is normal. Here 
the region J^ does not contain any straight lines and hence has 
at least one vertex. Indeed, if % lies in the plane (and, as noted in 
Section 2, such a case is possible), then .% is a convex polygonal 
region in the plane which contains no straight lines and, as explained 
in Subsection 2° of Section 5, must, therefore, have vertices. If the 
region A does not lie in the plane, however, consider its boundary. 
This consists of faces each of which, being a convex polygonal 
region containing no straight lines, must have vertices; and it is easily 
seen that a vertex of any face is simultaneously a vertex of the 
region #. 

Converging in the vertex A of the region # are at least three 
boundary planes for which the point A is the only common point. 
Indeed, if this were not the case, all the boundary planes passing 
throught A would either coincide or have a common straight line. 
But in that case a sufficiently small segment passing through A and 
lying in the common boundary plane or on the common boundary 
line would belong to .#, which is contrary to the definition of a 
vertex. 

This compels us to make some obvious changes in the method of 
finding vertices described in Subsection 2? of Section 5. That is, 
now we should say that a regular subsystem is a subsystem of 
three equations of the system 


ax + byy+c,z +d, =0 
WI RCN VS emere. (4) 
mX + Duy + Enz + dm 


I 
© 
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rather than two. provided the solution (x, y, z) of the subsystem is 
unique. When a regular system is understood like this, the method 
of finding vertices remains exactly the same as before, namely: 

To find all the vertices ofthe region X, one should find the 
solutions of all the regular-subsystems of the system (4) and pick 
out those which satisfy the original system (1). 

The theorem of Subsection 2° of Section 5 also remains valid: 
the changes to be made in the proof of the theorem are obvious. The 
remark that a normal system has no solutions if the region Jf 
has no vertices also remains valid. 

Example 1. Find the vertices of the region .# given by the system 
of inequalities 


2X tay) es LA 


x+2y+ z—101z0 (5) 
x+ y+2z-120 
ees 2 LO 


Here the corresponding homogeneous system of equations is 
of the form 


2x+ y+ 220 
X42y4 2=0 
x+ y+2z=0 
x+ y+ z-0 


Solving it we see that the unique solution is (0, 0, 0); so the 
system (5) is normal. 

To find the vertices we shall have to consider all the subsystems 
of three equations of the system (4): 


, 


BASE 'z—l-0 2x4 y+ z—120) 98 


x+2y+ z—120 X+2p+ 2—1=0 
Niet mse ei x+ y+ z—120 
2x+' y+ z—120 x+2y+ z2-1=0 
x+ y+2z-1=0 x+ y+2z-—1=0 
x+ y+ 2-1=0 ara Bt a) 


On making the required calculations we find that all the subsys- 
_fems are regular and that their solutions are the points 


(1/4,1/4.1/4), (0.0.1, (0, 1, 0), (1, 0, 0). 
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Que. 


of which the first does not satisfy the system (5) and the remaining 
three do. Therefore, the vertices of the region % are the following: 


A,(1, 0, 0), A2(0. 1, 0), A3(0, 0, 1) 


2°. The normal homogeneous system of inequalities (2). Each of the 
inequalities (2) determines a half-space whose boundary plane passes 
through the origin of coordinates. 

Here a single point, the origin of coordinates, is the intersection 
of the boundary planes (the system (2) is normal!). In other words, 
the set %o, the feasible region of the system (2), is a convex polyhed- 
ral cone with a single vertex. It follows from the enumeration of 
convex polyhedral cones given in Section 4 that in our case % 
is either an infinite convex pyramid, a flat angle, a ray or, finally, 
a single point (the origin of coordinates). Omitting the last case for 
the time being, we have in all the remaining ones 


H = (Bi, By, ~ B) 


where Bi, Ba, ... B, are some points chosen one at a time on 
each edge of the cone J'y (see Theorem 2 of Section 4). One can 
find such points from the following considerations. Each of. them 
(i) belongs to Jo, i.e. satisfies the system (2), and (ii) belongs to the 
line of intersection of two different faces, i.e. satisfies disproportion- 
ate* equations of the system (3). 

If it is found that the only point satisfying conditions (i) and 
(ii) is (0, 0, 0), then the region Jf coincides with the origin of 
coordinates. 

Example 2. Find the region .%9 of the solutions of the system 


2x y+ z20 
x+2y+ zz0 
x+ y+2z>0 (6) 


and, further, the region -4^ of the solutions of the system in Exam- 
ple |. , 

Notice first of all that the system (6) is connected with the system of 
inequalities (5) of Example 1; namely (6) is a homogeneous system 
corresponding to (5). Hence the system (6) is normal. 


+ Wecall the equations ax + by + cz = Qand a'x + b'y + cz = 0 "dispro- 
portionate" if at least one of the equalities a/a' = b/b' = c/c fails to hold: 


in this case the corresponding planes intersect along a straight line. 


47 


Here the system of two disproportionate equations can be set up 
in six different ways: 


x+2y+ z= | 2x yb z= j 2x yc - 


x+ y+2z=0 x+ y+2z= x+y+ z=0 
2x+ y+ z= x+2y+ z= Mo 
x+2y+ z= Dee yt z= x+y+ z=0 


For each of these six systems, choose two nonzero solutions: (x; y, v) 
and (— x, — y, — z). For example, one can take (3, — 1, — 1) and (— 3, 
1, 1) for the first system: inequalities (6) are satisfied only by the 
first of these solutions. Hence we have the point B, = (3, — 1, — 1). 
Proceeding in similar fashion with the remaining five systems we 
find the points B, = (— 1, 3, 1) and B, — (— 1, 1, 3). So the 
region Xo consists of points of the form 


UB, + t2Bz + (B3 = 
= (3t; — t2 — t3}, — t + 3t2 — 3,- t — fz + 3t3) 


where 1;, t5, t3 are arbitrary nonnegative numbers. 

We now turn to the system of inequalities (5) of Example 1: 
As already noted, it is (6) that is the corresponding homogeneous 
System. Therefore, the region J^ is of the form 


(Ai, As, A35 + Ao 
and consists of the points 
84i. 5242 $345 + GB, + GB) + (B, = 
=s; (L, 0, 0) -- s,(0, 1, 0) + s4(0, 0, 1) + 
+8, —1, —1)--t(—1, 3, = l)-t(—1,-1 39 
= (5-35 — tp = ts, 5) — ty + 315 = 13, $3 — tp — 15 JE 373) 


where the numbers li, f2, t3 are arbitrary nonnegative ones and 
St, $2, 53 are nonnegative and sum to unity. 

3°. The case where the System of inequalities (1) is not normal. This 
means that the feasible region # of the homogeneous system of 
equations (3) contains points different from the origin of coordinates. 
Since & is the intersection of planes, two cases are possible: 

l. Z^ is a straight line. According to Lemma | the region # 
together with each of its points P contains a straight line P + Z. 
Consider some plane 7 which is not parallel to Y. If we know 
which points of the plane 7 belong to the region ^ — denote 
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the set of these points by #;—we shall be able to find the 
region .# itself, for then KH=H>+L. 

But, whatever the straight line 2 may be, one can always choose one 
of the coordinate planes, XOy, xOz or yOz, as the plane ZF 
not parallel to it. Assume, for example, that 2 is not parallel to the 
yOz plane. Take this plane to be .Z. Then the set HK z — denote 
it now by #,.—is the part of the yOz plane contained in J^ 
(Fig. 44). To find this set, one should put x = 0 in the system (1). 
We then have the system of inequalities * 


Fig. 44 


biy+cız+d, >20 
bay + cuz + dy > 0 
which can be solved using the, method of Section 5. 
On finding the set .#,. we shall be able to write down 


X =H, +2 (8) 


(if the straight line . is not parallel to the yOz plane), 
which is the complete description of the region .%. 

Remark. If it is found that the set J,. is empty, so is X. 
This means that the system (1) is incompatible. 

Example 3. Find the region J^ of the solutions of the system 


—2x+ y+ 2-120 
—3x— y+4z-120 (9) 
— x—2y-43z 20 


* Jt is easily seen that the system (7) is now normal (see a similar 
footnote on page 42). 
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Consider the corresponding homogeneous system of equations 
tbe pte gae) 
—3x— y+4z=0 (10) 
— x—2y+3z=0 

Solving it we find that the third equation is a consequence of the 


first two, so that the system is reduced to the first two equations. 
The set of its solutions ¥ is the straight line along which the planes 


and 
—3x—y+4z=0 
intersect. 

Choose some point B on the straight line 2^ different from the 
origin of coordinates. To do this, it is enough to find some 
three numbers x, y. z (not vanishing all together) which satisfy 
the first two equations of the system (10). Take a d 1, for 
example. So £ is the straight line OB, where B = (1, 

It is easy to see that the straight line 2^ is not uen D. say. 


the yOz coordinate plane. Setting x — 0 in the system (9), we have 
the system 


yt z—1z0 
—-y+4z—1>0 
—2y + 3z >0 
in two unknowns y and z, which is normal. One can find its 


feasible region J/^,. using the method of Section 5. On making the 
required calculations we find that Jf. is a set consisting of a 


single point A(3/5, 2/5) (in the yOz plane). Hence the sought region 
HA consists of all the points of the form 


A+tB=(0,3, J+e(, 1, D= (534624 1) 


where t is any nonnegative number (the region J^ is a straight 
line parallel to 2). 

2. V? is a plane. Then take as a secant set 7 some straight line 
not parallel to this plane: in particular, one can take one of the 
coordinate axes. Assume, for example, that the z-axis is not parallel 
to 2^ take it to be 7. To find the set J^, the part of the 
z-axis contained in J^, one should put x = 0, y= “0 in the system (1 ) 
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We then have the system of inequalities 


cz 4d, >0 

mete ae ae i (11 

Cuz + d, = 0 ). 
which is solved without much difficulty *. On finding the set £^, 

we shall be able to write down (Fig. 45) 


cO el (12) 


(if the plane & is not parallel to the z-axis), which is the complete 
description of X. 


Remark. If the set X. is found to be empty, so is Jf. 
The system (1) is then incompatible. 7 
Example 4. Find the region .# of the solutions of the system 


T 


=—x+y—24+220 a5 


Here the corresponding homogeneous system of equations is 


of the form 
weha 
* 14 
—x-y- ze i (49) 


* The system (11) is normal. 
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ion i i f the first; therefore, the 
The second equation is a consequence o : : [ 
feasible region of the system (14) is the plane ¥ given by the equation 


x—ytz-0 


It is easy to see that this plane intersects the z-axis in a single 
point and hence is not parallel to the z-axis. Find the set X.. 
Setting x — 0, y=0 in the system (13), we have the system 


from which it follows that 
—1<z<2 (15) 


So X is the set X. + V consisting of the points of the form 
(0, 0, 2) +(x, y; —x + y) = (x, y z— x +y) where x and y are 
arbitrary and z satisfies the inequality (15). 

We conclude this section by formulating two theorems which gener- 
alize the last two theorems of Section 5 to the three-dimensional 
case. The only change we should make for this purpose in the 
formulations of the above-mentioned theorems of Section 5 consists 
in substituting the word “space” for “plane”. 

Theorem. Any (nonempty) convex polyhedral region in space can be 
represented as the sum 


KAi, Aa, n A,» + (By, Ba, os B) 
Theorem. Any set of the form: 
Ai, Ass s Ap? + (Bi, Ba, ., Bj) 
in space is either the whole of space or some convex polyhedral 
region in it. 
The proofs of both theorems follow almost word for word those 


of the corresponding theorems in the two-dimensional case. Their 
elaboration is left to the reader. 


7. Systems of Linear Inequalities 
in Any Number of Unknowns 


In the foregoing sections we concentrated upon systems of inequa- 
lities in two or three unknowns. This was dictated mainly by two 
circumstances. Firstly, these systems are simple to investigate and 
allow one to keep entirely within the framework of “school” 
mathematics. Secondly, (which is more important in the present 
case) solutions of such systems have visual geometrical meaning 
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(as points in the plane or in space). However, more common in 
applications (in linear programming problems, for example) are 
systems of inequalities involving n >3 unknowns. To pass them 
by in silence would greatly water down the presentation of the 
question, Therefore, we shall try to at least outline the situation 
for any n> 3. 

To visualize linear systems of inequalities in n unknowns we 
must turn to the so-called n-dimensional space. 

We shall begin with the definition of relevant concepts, con- 
fining ourselves only to the most essential. 

A point of n-dimensional space is defined by an ordered set of n 
numbers 


called the coordinates of the point. Such a definition is motivated 
by the fact, fundamental to analytic geometry, that a point in 
the plane is described by a pair of numbers and that in space 
by a triple of numbers. In what follows instead of saying “the 
point M has the coordinates x}, x5. .. x," we shall write M = 
(Xis Xa) ss Xp) or simply M (xj. X2» -a Xn). The point (0, 0, ..., 0) is 
called the origin of coordinates or just the origin. 

We shall first of all point out what is meant by a “segment” 
in n-dimensional space. According to Section 1, in ordinary space 
a segment M, M; can be described as a set of all points of the form 


syM, + 52M, 


where sı, s2 are any two nonnegative numbers whose sum is one. 
Going from three-dimensional space to n-dimensional space we adopt 
this description as the definition of a segment. More strictly, let 
Milx x2, 25 xDrande Mi" ct edi NE 
be two arbitrary points of n-dimensional space. Then we say that 
a segment M'M" is a set of all points of the form 
SM' + s"M" = 
= (s + S kis SX FSX, «4 SX, FISXA) (1) 


where s', s" are any two nonnegative numbers with a sum of 
one. When s'—1. s”=0 we have the point M'. when 
s —0, s"—1 we have the point M". These are the ends of 
the segment M'M". The remaining points (they are obtained when 
s 2 0, s” 7 0) are called the interior points of the segment. 


Of other concepts pertaining to n-dimensional space we shall need 
the concept of hyperplane. This is a generalization of the concept of 
plane in ordinary three-dimensional space. Here the prefix “hyper 
has quite definite meaning. The fact is that in n-dimensional 
space various types of “planes” are possible, viz. one-dimensional 
“planes” (they are called “straight lines”), two-dimensional “planes”. 
etc., and finally (n — 1)-dimensional "planes": it is the last that are 
called “hyperplanes”. 

Definition. In n-dimensional space. a set of points M (X1, X2 x,) 
whose coordinates satisfy the equation of the first degree 


MX, EOX +... +4,X, b =0 (2) 


n^n l 


where at least one of the numbers ay. a2... a, (coefficients of the | 
unknowns) is different from zero, is called a hyperplane. For n = 3 the 
equation (2) assumes the form a,x; + à3X5 + 43X3 +b = 0. which 
is none other than the equation of a plane in ordinary space 
(where the coordinates are denoted by x,. x5, X3 instead of x, y, Z 
as usual). 

With respect to the hyperplane (2) the whole of n-dimensional 
space is divided into two parts: the region where the inequality 


X, + aX +... +4,x, +b>0 (3) 
holds and the region where we have 
dX, Hax; Fax, +b<0 (4) 


These regions are called half-spaces. Thus every hyperplane divides 
the whole of space into two half-spaces for which it is a common 
part. 
The concept of a convex solid is also generalized to the n-dimen- 
- sional case. A set of points in n-dimensional space is said to be 
convex if together with any two of its points M' and M” it contains 
the entire segment M'M" as well. | 
It is easily shown that any half-space is a convex set. Indeed, 
let the points M’ (xz. x x) and M” (xf, x3... x") belong to the 
half-space (3). We prove that any point M of the segment M'M" 
also belongs to that half-space. 


The coordinates of a point M are written in the form (1) or, equiva- 
ently, as 


x= sy F= sx 
o; 5% - (I — s)xz 


X sssx cp (s 


Substituting into the left-hand side of (3) we have 
ai(sx + (1 — s)x7) + a(s + (1 — s)xz) + ... 
c asx, + (1 — sx”) + b= 
—os(aix + a2% +... + a,x’) + 
F (1 — s)(ayxf + asx2 +... + a,x!) + sb + (1 — s)b 
(we have replaced the number b by the sum sb +(1—s)b), 
which is equal to 


slax b ax +b] + — saxi + + a,x" + b] 


Each of the two sums in square brackets is nonnegative since 
both points M' and M" belong to the half-space (3). Therefore’ 
the entire expression is also nonnegative (for s > 0 and (1 — s) > 0). 
This proves that the point M belongs to the half-space (3), 
ie. that this half-space is convex. 

From the above it is easy to understand what geometric 
terminology should correspond to a system of linear inequalities 
in n unknowns. Given the system 


aX, + aX + t ax, az O 
byx, + baxa + + b,x, +b > 0 


(5) 


C4X1 Feka +... + Ox, +e 20 


each of the inequalities determines a certain half-space while all the 
inequalities together determine a certain region J^ in n-dimensional 
space which is the intersection of a finite number of half-spaces. 
The region Jf. is convex, since so is any of the half-spaces that 
form it. 

In n-dimensional space, by analogy with the three-dimensional case, 
a region which is the intersection of a finite number of half-spaces 
is called a convex polyhedral region or, if the intersection is a bounded 
set, simply a convex polyhedron. Here the words "bounded set" should 
be understood in the sense that the coordinates of all the points 
of a region under consideration do not exceed in absolute value 
a certain constant c: 

xil S c |x| & c for all points of a given region. 

"Thus, in n-dimensional space, a set of points whose coordinates satisfy 
the system (5) is the convex polyhedral region A” resulting from intersec- 
tion of all half-spaces corresponding to the inequalities of the given 
system. 
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Note again that if the region is bounded, it is called a convex 
polyhedron. A. A : 

The methods of actual description of the region J^ which have 
been considered for systems in two unknowns in Section 5 and for 
systems in three unknowns in Section 6 can, with appropriate changes. 
be extended to the case of n unknowns. We shall not go into 
this, however, for it would require rather much space. Besides, with 
a large number of unknowns these methods are not efficient, 
requiring as they do, too cumbersome calculations. 

It is remarkable that the general theorems on the structure of 
convex polyhedral sets in three-dimensional space remain fully valid 
in n-dimensional space, although requiring more elaborate proofs. 
Weshall restrict ourselves only to the formulations of the theorems and 
to the essential explanations. 

Theorem 1. The convex hull of any finite system of the points Ay, Aa. 
A, is a convex polyhedron. , 

For better intelligibility we shall emphasize that here a relation is 
established between two differently defined types of sets, viz. 
between the convex hull of a system of points Ai, As. A, which 
is designated as <A}, A5, ., A,> and defined as a set of all 
points of the form 


51344 $545 +.. + SA, 


where sy, Sams s, are any nonnegative numbers whose sum is one, 
and convex polyhedra, i.e. bounded regions resulting from inter- 
section of a finite number of half-spaces. 

In two- and three-dimensional spaces the validity of Theorem 
| is obvious (at least from the visualization of the convex hull). 
while in the multidimensional case it is not obvious at all and re- 
quires to be proved. 

Theorem 1' (the converse of Theorem 1). Any convex polyhedron 
coincides with the convex hull of some finite system of points. 

In fact we can state even more: a convex polyhedron coincides 
with the convex hull of its vertices. The definition of a vertex is the 
Same as in the two-dimensional case (a vertex is such a point of à 
polyhedron which is not an interior point for any of the segments 
belonging wholly to the polyhedron). It can be shown that the 
number of vertices is always finite. 

Theorem 2. Any set of the form (By, Bo,..., B,) either coincides 
with the whole of space or is some convex polyhedral cone with 
vertex at origin. i 

We remind that the symbol (B,, B, 


. B,) designates a set of all 
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points which can be represented in the form 
GB, + t3B3 E LB, 


where t4, £5... t, are any nonnegative numbers. A convex polyhedral 
cone is defined as the intersection of a finite number of half-spaces 
whose boundary hyperplanes have a common point (the vertex of the 
cone). The validity of Theorem 2 in three-dimensional space was 
established in Section 4 (Theorem 1 of Section 4). 

Theorem 2’. Any convex polyhedral cone with vertex at origin can 
be represented as (By, By,..., B,). 

The validity of the theorem for the three-dimensional case was 
proved in Section 4 (Theorem 2 of Section 4). 

Theorem 3. Any convex polyhedral region can be represented as 
the sum 


(Ay Ages Aj) + (By, Bana B) 


Theorem 3’, Any sum of the indicated form is either the whole 
of space or some convex polyhedral region in it. 


8. The Solution of a System 
of Linear Inequalities 
by Successive Reduction 
of the Number of Unknowns 


From his elementary algebra course reader is familiar with the 
method for solving systems of linear equations by successive reduction 
of the number of unknowns. For systems involving three unknowns 
x, y, z the essence of the method can be described as follows. 

From each equation of a given system, one finds the unknown z 
in terms of the unknowns x and y. The obtained expressions 
(containing only x and y) are then equated to one another. This 
leads to a new system, a system in two unknowns x and y. For 
example, the following was the original system 


2x —-3y z=-1 
x— yt2z-5 (1) 
4y- 2=5 
then, on solving each of the equations for z we have 
ze— 2x4 3y—]1] 
zx ly (t) 
z= 4y— 5 


after equating the expressions in the right-hand sides (it is enough 
to equate the first of the expressions to each of the rest in turn) 
we have the system 


1 1 5 
2x Sy — 1 aX *3Y 5 o 
—2x+3y-1= 4y—5 


in two unknowns x and y. The system (2) can be handled in the 
same way. From the system of equations (2) we find 


Hb e. 7 
Ti pni 


y= —2x4+4 


after which we get the equation 
See EAS 
sx tas —2x+4 


involving now one unknown x. On solving the equation we find 
x= l, after which from the system (2) we find y = 2, and finally 
from the system (I) we have z= 3. É 

It turns out that something of the kind can be done with any 
system of linear inequalities. Solving a system of inequalities in 7t 
unknowns x. x, ,, x, is thus reduced to solving a system of 
inequalities in n — 1 unknowns x;.... X, p then the resulting system 
can be reduced to a system in n — 2 unknowns x, x, and so 


An 
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Consider any of the inequalities of the system (S). It is of the form 
AXi 3X3 Hae Han daz (3) 
If a, = 0, leave the inequality unchanged. If a, <0, transfer the 
term a,x, to the right-hand side of the inequality and divide both 
sides by a positive number — a, to get the inequality of the form 
bixi +. + bai Xn bx, 
In the case of a, > 0. transfer all the summands but a,x, to the 
right-hand side of the inequality and divide both sides by a, to 
get the inequality 
Mee CL eo Chen en leat 
So on multiplying each of the inequalities of the original system 
by an appropriate positive number we get the equivalent system 
of the form 


(T) 


where ^P, Py Qi Op Ris Re áre expressions of the form 
jus Ps Q, 


nint MY 
diXj bet daan- + 0 (not containing x,)*. 


*Of course, if the system (S) has no inequality in which a, < 0, then the 
system (T) will not have the first block. Similarly, if there are no inequali- 
ties with a, 7 0. then the sy: (T) will not have the second block. - 
Finally. if (S) has no inequalities with a, = 0, then the third block will 


be missing from (7). 7 
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In shorthand notation the system (T) can be written down as 


P, > xX, > s 
R,>0 
where « is any of the numbers 1, 2... DP. B is any of the numbers 
1, 2... q, and y is any of the numbers 1, 2.., r, 
Along with the system (T) consider the System 


P, > Qs (S) 
R,z0 
(where o is any of the numbers 1. 2- D. is any of the 
numbers 1, 2.. q, and y is any of the numbers 1, 2,., 7) 
in n— 1 unknowns Xise Xq2.*. Let us agree to call this 
system concomitant (with respect to the original system (S) 
or to the equivalent system (T)). 

There is a close connection between the solutions of the systems 
(S) and (S’) which is expressed in the following 

Theorem. If the value of the last unknown X, is discarded from 
any solution of the system (S), then we get a certain solution of the 
concomitant system (S'). 

Conversely, for any solution of the concomitant System (S') it is 
possible to find such a value of the unknown x, on adjoining which 
we get a solution of the original system (S). 

The first statement of the theorem is obvious (if a set of values 
of the unknowns satisfies the system (S), then it does the system (T) too; 


but then all the inequalities of the system (S") hold for the same set too). 
We prove the second statement. 
Let 


DDA " BG 
Ej S Xe DUI ES SUM 


be a solution of the system (S'). Substituting into the expressions 
Pis Ps Orban Qp Rise R. we get some numbers p... 
P MONERON RD R?. The following inequalities must hold 
for them ` 

„P> OF (4) 
(x being any of the numbers Lh. 2:55 p. and f any of the 


* Ifthe first or the second block is found to be missing from the system (7). 
then (S') will consist only of the inequalities R, >0. If the third block is 
missing from (T), then (S’) will have only inequalities PL > Qy 
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numbers 1, 2.., q) and 
R$ 2z0 > (5) 


(y being any of the numbers 1, 2...., r). 

The first group of the above inequalities (i.e. (4)) shows that each 
of the numbers Q°.... 0 is not greater than any of the numbers 
PY. P5. But in such a case there must be a number xo which 
lies between all the numbers Qj... Q? and all the numbers P9, 
aes 


p> 


PS > xa > O8 


where œ is any of the numbers 1, 2...., P, and D is any of the 
numbers 1, 2.., q (Fig. 46). These inequalities together with the 


Th 
— H 
a p o s o 0 po 
$54 7H RE 
Fig. 46 


inequalities (5) mean that the set of the values of the unknowns 
2 E) n E) 4 4 
Xi S Xis Xn-1 = Xn 1, XQ XD 


is a solution of the system (T) and hence of (S). Thus the theorem 
is proved. 

The following two additions to the theorem will play an important 
role below. 

1. A system (S) of linear inequalities is compatible if and only if so 
is the concomitant system (S). This is a direct consequence of the 
theorem proved. 

2. All the solutions of the original system (S) can be obtained in 
the following way. One must adjoin to each solution x9, .., xo Of, 
the concomitant system (S') any of the numbers x? lying between all 
the numbers Q3, .., Q9 and all the numbers P9, ..., P). As a matter 
of fact this proposition was proved during the proof of the theorem. 

We shall say a few words about the geometric meaning of the 
theorem proved above. Suppose that (S) is a system of inequalities in 
three unknowns x, Y, z. The concomitant system (S’) is a system 
in two unknowns x, y. Denote by X (S) the feasible region of 
the system (S) (it is a certain set of points in space) and by # 
(S") the feasible region of the System (S") (a set of points in the plane). 
In geometric terms the theorem proved has the following meaning. 

The region ^ (S) is the projection of the region X (S) on the 
XOy coordinate plane. 
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So for an arbitrary system (S) of linear inequalities involving 
unknowns x4. XA, we have constructed a new, concomitant system 
(S) involving the unknowns -x;. Xa, .. x, ,. and for the system 


(S‘) one can in turn construct the concomitant system (S") (in- 


volving the u. Owns X1, X5. .... X4. 5). for the latter system one can 
construct the concomitant system (S"') and so on. 
Continuing the process we shall. after a number of steps, come 
to the system (S" - ") consisting of inequalities in one unknown X 
It follows from Proposition | stated above that the system (S) is 
compatible if and only if so is the system (S"- '), and answering 
the question concerning the compatibility or otherwise of a system 
in one unknown presents no difficulty. Thus it becomes" possible 
for us to learn by quite simple calculations whether the system (S) is 
compatible or not. 

Suppose that the system is compatible. Then the problem arises 
of solving the system or, speaking at greater length, of indicating 
all of its solutions (all sets of values of the unknowns which satisfy 
the inequalities of the given system). We shall adopt the point of 
- view (though it may at first appear to the reader strange) that the 
system (S) is solved if the systems (S^) (S"). .... (S"7 ) are constructed. 
We shall presently see what explains such a point of view but first 
we shall introduce a 

Definition. A set of values of the first k unknowns 


: Xie e. xf 


is .called feasible if it can be extended to the solution of the 
original system (S), i.e. if there exist numbers x2... .... x? such that 
the Set xf, Ln xp. x24, 0. x9 ds the Solution of the system (S). 
Once the systems (S^), (S") and so on are constructed, we are in 
à position: 

(s iy find all feasible values of the 


2) to find for any particular feasible value xt all compatible 
values of the unknown X», such that together with x? they form a 
feasible set (they are found by substituting X? into the system (s i) 

3) to find for any particular feasible set Xs x2 call compatible 


values of the unknown X3 (they are fou ituti ? 
) 3 nd by subs gx 
and x into the system ($" 73)) and so on. i RT 


It is in this sense that one should understand our statement 
that the system (S) is solved if the systems (S'), (S") (Se) 
are constructed. NE wa 


unknown x, (from the system 


[ 
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Example. We solve in the above sense the system 
Dao O. 


On solving each of the inequalities for z we write down the 
system in the form 


The concomitant system is of the form 
jx-y-22 2x—-3y4]1 
jx-y-22 —5x4 y—2 
—x—y+t4> 2x—3y-cl 
—x—y+42—-—5x+ y—2 
or, after the cancelling of the terms, 
x+4y—320 
>0 


—3x42y 4320 
4x—2y +620 


On solving each inequality for y we write down the system in 
the form 


3 
as 
24 


ak 


vv WV 


Y 
a 
# 
x w 
SAS 
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The concomitant system is of the form 


2x32 —$x4i 


X32 3x-3 


xz0 
it is equivalent to one inequality 
x20 (8) 


Thus the original system is compatible. According to the point of 
view we have adopted the systems (8), (7), (6) provide the solution 
of the set problem; namely, inequality (8) shows that there exists a 
solution (x, y, z) of the original system involving any nonnegative 
X. If a particular value of x is chosen, then from the system (7) 
one can find the feasible values for y. If particular values for 
x and y are chosen, then from the system (6) it is possible to 
find the feasible values of z. 

Set x = 1, for example; then from the System (7) we have the 
following inequalities limiting y: 

52y2i 
Take, for example, y = 4. Setting x= 1, y = 4 in the system (6), we 
have the following inequalities limiting z: 


yz 22 
zu i 
z2—9 
zz —3 
or simply 
Sal >> — 3 
Setting, for example, z= — 2, we have one of the solutions of 
the original system: x = 1, y=4 z= 


9. Incompatible Systems 


i So far we have been dealing mainly with such systems of inequali- 
ties which have at least one solution (which are compatible). As 
far as incompatible systems are concerned, studying them may on 
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the face of it appear a needless pursuit; besides it seems implausible 
that one should be able to relate any interesting theory to such systems. 
In fact everything is quite different. The properties of incompatible 
systems are not only of interest in themselves, but also give key 
to the understanding of a variety of important facts. Thus the main 
linear programming theorem (the duality theorem, see Section 14) is 
in the long run derived from some properties of incompatible systems. 
Consider an arbitrary system of linear inequalities. For notational 
convenience consider for the time being that the number of unknowns 
` is three, although everything stated above equally applies to systems 
in any number of unknowns. 
So we are given the system 


axtbytez+d,20 | 
z0 


aax + bay +22 + dz 
(1) 


aux + buy cuz + d, > 0 
Multiply both sides of the inequality of (1) by a nonnegative 
number K;, both sides of the second inequality by a nonnegative 
number k, and so on, and then add the obtained inequalities 
together to come as a result to the inequality 
(kiai + kaa +... + k,a,)x + 


+ (Gb, + kaba +. + kub,)y + 
+ (kici + kaca + + kmen) + 
+ kydy + kada +...+k,d,, Z0 (2) 


which will be called a combination of the inequalities of (1). 
It may happen that a certain combination of the inequalities of 
(1) is an inequality of the form 


Oxx+0xy+0xz+d30 (3) 


where d is a strictly negative number (division by |d| then leads 
to the inequality — 1 > 0). It is clear that no set of values of the 
unknowns satisfies such an inequality, therefore in the case under 
discussion the system (1) is incompatible (it has no solutions). It is 
quite remarkable that the converse proposition is also valid; namely, 
if the system (1) is incompatible, then a certain combination of its 
inequalities is of the form (3) where d — 0. 

We shall now prove this proposition in the general form (i.e. for 
systems in any number of unknowns), but first let us introduce the 
following definition. Let us agree to say that the inequality 


ax + by+cz+d>0 
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is inconsistent if it is not satisfied by any set of values of the unknowns. 
Obviously, any inconsistent inequality is of the form (3) where d < 0 
(prove this!). The proposition to be proved can now be formulated 


as the following theorem. à N 
Theorem on incompatible systems of inequalities. If a system of 


linear inequalities is incompatible, then a certain combination of these 
inequalities is an inconsistent inequality. 2 

" The proof is made by induction on n, the number of unknowns 
in our system. 

When n=1 the system is of the form 


a,x+b,>0 


(4) 


a,x + b,, > 0 


One can suppose that all the coefficients a, a5,.., a,, are nonzero. 
Indeed, if for example, a, = 0, then the first inequality is of the form 
0x x b, Z0; if the number b, is nonnegative, then such an ine- 
quality may be discarded; if b, is negative, however, then the first 
inequality of our system is inconsistent and there is nothing to prove. 

So let us suppose that none of the numbers 04, d5,.. Ay, IS ZETO. 
It is easy to see that these numbers are sure to contain both positive 
and negative numbers. Indeed, if all of the indicated numbers had the 


same sign, were they positive, for example, then the system (4) would 
be reduced to the form 


b, 
x2z—— 
ar 
bz 
x2-— 
a 
bm 
xz—-— 
a 


and would hence be compatible. 
Assume for definiteness that the first k numbers of 0,, d5,..., Ay ATE 
positive and the remaining m — k are negative. Then the system (4) 
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O T cA uo o ————— ÓÉ——— 


is equivalent to the system 


(5) 


Choose the largest among the numbers — (b,/aj),.., — (b,/a,); let 
it be —(b,/a,), for example. Then the first k inequalities of the system : 
(5) can be replaced by a single (first) inequality. Similarly choose 
the smallest among the numbers — (bkg ./a,. Doon — (b,/a,), the 
remaining m — k inequalities of the system (5) can then be replaced 
by a single (last) inequality. Thus the System (4) is equivalent to 
the system of two inequalities: 


zb 
ay 
bm 


Am 


xz 


and its incompatibility means that 
b, b 
-—>-— (6) 
From (6) we have å 
b,a, — bya, < 0 (7) 


(it should be remembered that a, > 0 and a,, < 0). If the first inequal- 
ity of (4) is now multiplied by a positive number — a, and 


the last inequality by the Positive number a,, then after addition 
we have the inequality 


0x xc (bna; ae bian) > 0 (8) 
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which is by virtue of (7) inconsistent. For systems in one — 
the theorem is thus valid. We now assume that the statement of the 
theorem is valid for systems in n— 1 unknowns and under this 
assumption establish its validity for the case of n unknowns. 

So suppose an incompatible system of linear inequalities in n un- 
knowns X;, X2, .. x, is given; keeping to the notation of the previous 
section we call it "system (S)". We construct the concomitant system 
(S’) for it; the latter will be incompatible after (S). Since the:number 
of unknowns in the system (S’) is n — 1, the hypothesis of induction 
applies to it. This means that a certain combination of inequalities 
of the system (S’) is an inconsistent inequality. But it is easily seen that 
each inequality of the system (S') is a combination of inequalities 
of (S): indeed, if we simply add up the inequalities P,2X, and 
x, > Qg of the system (S), we get P, + x, > x, + Qp or P, > Qp i.e. one 
of the inequalities of the system (S’) (see Section 8, systems (T) and 
(S’) ). Hence it follows in turn that a certain combination of inequali- 
ties of the original system (S) is an inconsistent inequality. Thus 
the theorem is proved. 

The theorem on incompatible systems of linear inequalities is 
just one of the manifestations of the deep analogy existing between 
the properties of systems of linear inequalities and those of systems 
of linear equations. Let us try to substitute the word “equation” 
for “inequality” in the formulation of the theorem. We shall have the 
following proposition: 

Ifa system of linear equations is incompatible, then a certain combi- 
nation of these equations is an inconsistent equation. 

It turns out that this proposition is also valid. It is called the 
Kronecker-Capelli theorem in a somewhat different formulation and 
is proved in linear algebra (this is the name of the part of mathema- 
tics which studies linear operations, i.e. such operations as addition of 
points or multiplication of a point by a number in n-dimensional 
space). For the above to be better understood, however, it is necessary 
to introduce clarity into the concept of combination. A combina- 
tion of equations is constructed in the same way as a combination 
of inequalities, the only difference being that one is allowed to 
multiply the given equations by any, not only nonnegative, numbers. 
As in the case of inequalities the word “inconsistent” refers to an 
equation having no solutions. It is easily shown that an inconsistent 
equation must certainly be of the form 


Oxx,+0xx,+..40xx,+5=0 (9) 


where b is a nonzero number (division 


of both sides by b gives 
the “equation” 1 = 0). LACE: 


68 


Especially important is one particular case of the theorem on in- 
compatible systems of inequalities, namely, the case where a given 
system contains the inequalities 


m0, a, 20, 24 x20 (10) 


Denoting the rest of the system by (S), one can say that the problem 
is to find all nonnegative solutions of the system (S) (i.e. the solu- 
tions satisfying conditions (10)). If the problem has no solutions, then 
according to the theorem just proved a certain combination of in- 
equalities of the system (S), 


aX, +X. + aX, +a>0 (11) 


yields in conjunction with a certain combination of inequalities of 
(10), 


kyxy + kaxa +... +k,x, 20 
(kis k2,.., k, are nonnegative) 
the inconsistent inequality 
Oxx+0xx.+..4+0xx,+c20 
where c is a nonnegative number. Hence 


a, = —k, <0, ag = — k3 < 0, ..., a, = —k, <0, 
a<0 


We shall formulate the result obtained as a separate proposition. 

Corollary of the theorem on incompatible systems. If a system of 
inequalities has no nonnegative solutions, then a certain combination 
of these inequalities is an inequality of the form (11) where all the 
coefficients ay, az~, A,<0 and the absolute term a « 0. 


10. A Homogeneous System 
of Linear Inequalities. 
The Fundamental Set of Solutions 


In Section 8 we have discussed a method for finding solutions 
of systems of linear inequalities. In spite of its having many obvious 
merits the method fails to give answers to some questions: for 
example, it does not allow one to review the set of all solutions 
of a given system of inequalities. It is to this end that this and 
the next section of our book are devoted. The main difficulties, it 
will be seen, arise in considering homogeneous systems which are dealt 
with in the present section: the general case (ie. the case of a 


69 


‘nonhomogeneous system of inequalities) is examined in Section 11. 
There is no need here to confine ourselves to the case of two or 
three unknowns, from the outset we shall consider a system of any 
number of inequalities in any number of unknowns. For the reader's - 
convenience our account is broken down into a number of points. 

1°. A linear function of n arguments. The general form of a ho- 
mogeneous inequality in n unknowns is 


1X; + d3X; +... ax,20 
Consider separately the expression in the left-hand side of the in- 
equality 
01X4 + 3X3 t + a,X, (1) 
It is called a linear function. The role of the arguments is played by 
n variables xi, x3,..., x,. One can assume, however, that function 


(1) depends on one rather than n arguments; this argument is 
the point 


X = (Xx, Xa. x) 


of n-dimensional space. 


Let us agree from now on to designate function (1) as L(X) for 
short: 
< L(X) = ax; + aoxa +... + a,x 


n^n 


if we are given several functions like this instead of one, we shall 
designate them as L(X), L(X), etc. 


Nn establish the following two properties of a linear function. 


L(kX) = kL(X) 


where X is any point and k is any number. 
2 


L(X + Y) = L(X) + L(Y) 
where X and Y are any two points. 
Property 1 follows in the obvious way from the equality 


ay(kx4) + az(kx2) + ... +.a,(kx,) = 


= K(ayx, + aaxo +... + a,X,) 
We now prove property 2. Let 


X= (x1, X2 vy x,) and Y= (i. Y2; +5 y,) 
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Then . 
L(X + Y) = a(x; + yi) + a(x + y3) + 
cr A(X, Vy) = (aixi + daxa + ~- + AX) + 
+ (ayy, + a22 +... + a,y,) = L(X) + L(Y) 


2°. Some properties of the solutions of a homogeneous system 
of linear inequalities. Given a homogeneous system of m linear in- 
equalities in n unknowns 


dXi + daXa c +4,X, 20 (the first inequality) 
bx, + b2X2 +. + bX, > Q (the second inequality) 
CX, + C2X2 +. HCX 20 (the mth inequality) 


Designate the left-hand sides of the inequalities as L,(X), L;(X )..... 
L,(X) to rewrite the system as 


where 
X = (Xi Xz Xn 


We establish the following two propositions. 

1. If the point X is a solution of the system (2) and k is any 
nonnegative number, then the point kX is again a solution of the 
system (2). 

2. If X and Y are two solutions of the system (2), then X + Y is 
again a solution of the system (2). 

Both propositions follow readily from the properties of a linear 
function proved in sect. 1°. Indeed, suppose i is any of the numbers 


1, 2.., m. We have 
L(kX) = kL(X)2 0 
(for k z 0 and L(X) z 0) as well as 
L(X + Y) = L(X) - L(Y) 2 0 
(for L(X) 2 0 and L(Y) > 0). 


Propositions 1 and 2 immediately yield the following consequence. 
If certain points Xi, Xa... X, are solutions of the system (2). 
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then so is any point of the form 
kX, tX, +... +k,X, (3) 


where ky, ka... k, are nonnegative numbers. " 

Indeed, by virtue of Proposition 1 each of the points kX, 
k2X>,.., k,X, is a solution of the system (2) but then by virtue 
of Proposition 2 so is the sum of the points. Let us agree to call any 
point of the form (3), where k,, k;..., k, are nonnegative numbers, a 
nonnegative combination of the points x X3,.., X,. Then the above 
consequence will allow the following statement. 

A nonnegative combination of any number of solutions of the homo- 
geneous system (2) is again a solution of this system. 


3°. The fundamental set of solutions. We introduce the following 
definition. 


A set of a finite number of solutions 
XR Lx 


of the homogeneous system (2) is said to be the fundamental set 


of solutions if any solution X of the System can be given by the 
formula 


X-hX, tkhX; 4.24 kX, (4) 


where ky, kas., k, are nonnegative numbers. It follows that in this 
case formula (4) in which ki, ka, k, are any nonnegative numbers 
gives a review of all solutions of the system (2). Hence it is clear 
that the problem of finding the fundamental set of solutions is 
one of primary importance in the investigation of the system (2). 
The development of an algorithm which would allow the fundamental 
set of solutions for any system (2) to be found using quite simple 
operations is what we set ourselves as the final object. 

4°. Construction of the fundamental set for a sy. 
one inequality. Construct the fundamental set of 
inequality 


stem consisting of 
solutions for the 


4X1 + dx; +... + a,x, > 0 (5) 
where the numbers ay, a,..., a 


, do not vanish together, 
To do this, consider along 


with inequality (5) the equation 

4,X; + aX. +... + a,x, = 0 (6) 
The properties of a linear function Proved in sect. 1° 
the following two propositions: 


1. If X is some solution of equation (6) and k is any nuinber, 
then kX is again a solution of equation (6). 


readily yield 
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2. If X and Y are two solutions of equation (6), then X +Y 
is again a solution of equation (6). 

The proof of the propositions is left to the reader. 

According to the premises there are some nonzero numbers among 


ai, Ary, a, Set a, #0, for example. 
Then the equation can be written down as 


1 
x= a, e +azxa + + ap aXy2j) (6) 


Setting x = 1, x, =0...., x, , =0, we find from the last equation 
that x, = —(a,/a,). Thus the point 


a -(: Oey = a) 
ay 


is a solution of equation (6), Proceeding in a similar way one can 
obtain the solutions 


Now let 
X = (04, Oa, .., &4-1. On) (7) 


be any solution of equation (6). According to (6') we have 


1 
a = —— (ait + d385 +. +4, 0, 4) = 


a, 


= ENT zs e = Ge 
Considering the point 


Oy Xy + 2X2 + F0, Xn-1 = 
=a, (1, i. -2 ex (0. m -#)+. 


convinces us that its coordinates coincide with those of the point X. 
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Therefore the following equality is valid 


X = 0X, +X +... +O, 4X, (8) 
We now add one more point 
Xn = —UG + Xi us Xu) (9) 


to the points X;, X5. .., X, , constructed earlier. It follows from 
the properties of the solutions of equation (6) referred to at the 
beginning of this subsection that the point X, is also a solution. 
Now it is easy to prove the fact that any solution X of equation (6) 
is a nonnegative combination of the solutions Kay Xe; s X, Pete 
Indeed, let x be a positive number exceeding any of the 
numbers |o]. el. ds xnl- It follows from (8) and (9) that 


X — (a, - o) X, + (a. +a) X4 +... 
(x: -Fo)X,-, + «X, 
which proves our statement. 
For brevity of further writing denote the left-hand side of 
equation (6) by L(X). Choose some solution of the equation 
L(X) — 1 and denote it by X,«,. We state that the set of points 


BEd og Yep is Kies (10) 


n-\ 


is the fundamental set of solutions for inequality (5). 
Indeed, each of the points satisfies inequality (5). Now let X' 


be any solution of the inequality; hence L(X’)=a where a> 0. 
Then the point 


X= Xx AX ns} 
satisfies inequality (6) for 
L(X) = L(X) - aL (Xy41)=a—ax t=0 


If we now write down 
X'— X-raX,4 


and recall that the point X is a nonnegative combination of the 
points Xis. Xp- p X, it will be clear that X’ can be represented 
as a nonnegative combination of the points (10). 

Consider a specific example. Suppose it is required to construct 
the fundamental set of solutions for the inequality 


— 2x, — 4x. +4350 (11) 
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in three unknowns xj, x, xs. 
First of all we write down the equation 


— 2x, — 4x, + x3=0 
and solve it for one of the unknowns, X3, for example: 
X3 = 2x, + 4x; , 
We now successively set one of the unknowns X;, X3 (contained 


in the right-hand side of the equation) equal to 1 and the remaining 
unknowns to zero, and get the solutions 


X,=(1, 0, 2), X,=(0, 1, 4) 
As X; we then take the point 
X3= —(X,+X2)=(—1, —1, — 6) 
and finally take as X, one of the solutions of the equation 
— 2x, — 4x, +.x3=1 
X4 = (0, 0, 1), for example. 


The points X,, X2, X4, X, form the fundamental set of solutions 
for inequality (11). The general solution is of the form 


X = kX, +k2X2 + k3X3 + k4X4 = 
= ky(1,0,2) + k2(0, 1,4) + k3(— 1, — 1, — 6) + k4(0,0, 1) 
or ] 
xX, =k, — k 
X2 = kj — ks 
X3 = 2k, + 4k; — 6k3 + k4 
where k,, ka, ks, ką are arbitrary nonnegative numbers. 
5°. The rearrangement of the structure of the fundamental set of 


Solutions when another inequality is added to the system. In order to 
learn how to construct fundamental sets of solutions, first consider 


a problem like this. 
Suppose we are given the homogeneous system 
L(X)20 
L{X)>0 (2) 


L(X)20 


of linear inequalities. Suppose further that we know the fundamental 
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set of solutions 
May X25 MX 


set 0j 
for that system. It is required to construct the fundamental set of 


s iti her 
solutions for the system resulting from the addition of anotl 
inequality 


p 


L(X)20 (13) 
to (12). y 

The solutions of the system (12) are precisely all nonnegative 
combinations of the points: X,, Xa- X, We must pick among 
these combinations those which would satisfy inequality (13) and, 
what is more, constitute the fundamental set of solutions for the 
system (12), (13). " 

All points X,, Xa ..., X, can be divided into three groups 
with respect to the function L(X), the left-hand side of inequality (13): 
the points for which L(X) > 0, the points for which L(X) « 0, and 
finally the points for which L(X)=0. We designate the points of 
the first group as Xj, X$., X*. the points of the second 


group as Xi; X;.., X^, and those of the third group às 
AUR IE E 
Thus 


Eua Xf X Ki MP. x9 


s 

are the same points X,, X,, 
perhaps in a different order. 
All the points X7 (a= 1 


^» X, except that they are located 


and so do X?(y =1,..., s). “As to the points Xj (B = 1...., I) none 

of them is a solution of inequality (13). However, of each pair 
Xe, Xp 

(one “plus” and one “minus” 


point) one can form a nonnegative 
combination 


aX! + bX; (14) 


so that it should satisfy the condition L(X) — 0. To do this one 


should take 
a— —L(Xy) b= L(X7) (15) 
Indeed, the numbers a and b are positive and besides 
L(aX7 + bX) = al(X) + bL(X;,) = 
= L(X5)L(X;) T L(X7)L(X5) -0 
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Denote the point (14), where a and b have the values indicated 
above, by X5, 


X= — L(Xg)XZ  I(X7) Xy (16) 


The answer to the set problem is given by the following 
Theorem. The points 


1 Sine Aian s recs XD ONO o XE (17) 


(there are k +s+ kl points here in all) form the fundamental. set 
of solutions for the system (12), (13). 
To prove the theorem, we first establish the following lemma. 
Lemma. Any nonnegative combination of the points X? and Xy 
can be represented as a nonnegative combination of the points X;, 
X% or else as a nonnegative combination of the points X p Xap. 
Proof of the lemma. Let 


X —cX; t dXg 


be a nonnegative combination of the points X7 and X;. Along 
with X consider the point 


XAK bXy 
where the numbers a and b are given by formulas (15). Compare 
the two ratios, c/a and d/b. If the first is greater than the second, 
then, setting d/h = à, c/a =A + where tp > 0, we have 
X = (ka ua) X7 c AbXg =X% + paxt 


i.e. the point X can be represented in the form of a nonnegative 
combination of X7 and X$2, If the above ratios are equal, 
then X = AX2,. Finally, if c/a < d/b, then, setting c/a = à, d/b = 
4+ where p 0, we have 
X =X% + bX; 

Thus the lemma is proved. , 

Proof of the theorem. Notice first of all that each of the points 
(17) satisfies the system (12), (13). To prove the theorem, all that 
remains therefore is to check the following: if some point X is 
à solution of the system (12), (13), then it can be represented in 
the form of a nonnegative combination of the points (17). 

Being a solution of the system (12), the point X can be represented 
as a nonnegative combination of its fundamental points X 1, X 2... X. 

X —aXj-.caXé-qTbXp.- 
ww BX, Me X0 us ocX? (18) 


TI 


where all the coefficients a... aj, bis- bp Cis- c, are Honda d 

If all the numbers 5b,,., b, are zero, then there is obviously 
nothing to prove. Suppose therefore that there are strictly positive 
numbers among those indicated above. Note that among 4... Ay 
there are also strictly positive numbers, for otherwise we should have 


L(X)= bı L(X7) +.. +bıL(XT) + c,L(X9) +... 
ache X2) 


which is impossible since X satisfies inequality (13). m 

Suppose for definiteness that a; > 0 and b, > 0. Making use o be 
lemma we can replace the sum a, X1 +b,X7 by a Donner d 
combination of the points X}, X9, or of the points Xi. Xie » 
such a replacement is made in the expression for the point ? , 
then the total number of nonzero coefficients of X7... Xi 4 Xi ves 
X; will decrease at least by 1. If at the same time dts 
found that not all of the coefficients of X; e. Xr are zero in the 
newly obtained expression for X, then we again replace one of the 
sums of the form aX -- bXg by p EEN A combination of 
the points X4, X% or that of the points X5, X%; as a result the 
number of nonzero coefficients of Xý, X, Xj.., XT again 
decreases at least by 1. We continue in this way till an expression 
is obtained for the point X in which all the coefficients of Xj. 
Xy are zero. We then come to the equality of the form 


X-—aXi-c-.-aXi + Y dX, F 
ap 
tox) uoa cX9 


where all the coefficients on the right are greater than or equal to 
zero. But this is exactly the required representation for X. Thus 
the theorem is proved. 

6°. The existence and the method for the construction of the funda- 
mental set of solutions. Let us consider an arbitrary system of homo- 
geneous linear inequalities, For the first inequality of the system we 
can (using the method described in sect. 4°) construct the fundamental 
set of solutions. On adjoining to the first the second inequality 
we can, on the basis of the theorem Of sett. 5° construct the 
fundamental set of solutions for the system consisting of the first 
two inequalities. We then adjoin the third inequality and so on until 
we have the fundamental set of solutions for the whole of the original 
system of inequalities. This Proves the existence, and at the same time 


points out the method for the Construction, of the fundamental set 
of solutions. 
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1 


Of course, if in a given system of inequalities there is a subsystem 
for which one can at once point out the fundamental set of so- 
lutions, then that subsystem should be taken as the starting point; 
successively adjoining to it the remaining inequalities one will arrive 
after a number of steps at the desired fundamental set. 

Example. It is required to find for the system 


Li (X) = — 3x, — 4x2 + 5x4 — 6x4 20 
L,(X)= 2x, + 3x2 — 3x3 + x4 20 


all nonnegative solutions, ie. all solutions which satisfy the con- 
ditions 


(19) 


x,20 
x2>0 
x320 (20) 
x420 
To put it another way, it is required to find a general solution of 
the system (19), (20). 
It is easily seen that for the system (20) the fundamental set of 
solutions will be the set of the points 
X,=(1, 0, 0, 0), -X = (©, 1, 0, 0) 
X, = (0, 0, 1, 0, X4=(0, 0, 0, 1) 
(indeed, any solution (%1, %2, %3, %4) of the system (20) can be repre- 
sented in the form à, X, + 42X2 + 3X3 + o4 X 4). Adjoin the first ine- 
quality of (19) to the system (20) and, using the theorem of sect. 5°, 


construct the fundamental set of solutions for the system obtained in 
this way. For computational convenience make the following table: 


X, 1 0 0 0 5 
X, 0 1 0 0 -4 
X, 0 0 1 0 5 
X; 0 0 0 1 —6 
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i i row of the table is one of the fundamental 
"m atic (20) and the value of the function L, (x ) for 
that point. It is seen from the table that the only point of X; pe 
is X4. that the points of Xp type are X,, X» X4; and that 
there are no points of X9 type in this case. 
Find the points of X2, type. They are 


X9 —53X4 + 5X, = (5, 0, 3, 0) 
X52 =4X; + 5X, = (0, 5, 4, 0) 
X$,= 6X4 -- 5X4 = (0, 0, 6, 5) 
To avoid making further writing more complicated, designate these 
points as Y;, Y», Y, (respectively) and write Y3 instead of X3. 
The points Y, Y;, Y», Y; form the fundamental set of solutions 
for the system consisting of (20) and the first inequality of (19). 
Adjoin the second inequality of (19) to that system and make 
the following table: 


ot 
e 
o 


MS 
n 
2 


Pod 
o 
[i 


Faz 
o 
o 


It is seen from the table that the role of the points Y? is now 
played by Yi. Y, and that of the points Yy by Y, Y, and that 
there are no points of Y? type. 

We find points of Yop type: 


Y= 3X-c Y,-(15 0, 10, 0)= 


)75(3, 0, 2, 0) 
Y= 3% +3%=(0, 15, 15, 0)=5(0, 3, 3, 0) 
= BY + Y= (65, 045 5)=5(13, 0, 9 1) 
yh = 13¥, + 3Y, =( 0, 65, 70, 15)=5( 0; 18: 14, 3) 


The points Y,, Y2, Y. Y9, 
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» Yt. Y9, form the fundamental set 


ofsolutions for the system (19), (20). The general solution is of the form 
X —aY, + bY, + c¥i3 + dYs + ea +f Ya = 
— a(5, 0, 3, 0) - b(0, 5, 4, 0) + 5c(3, 0, 2, 0) + 
+ 5d (0, 3, 3, 0) + 5e(13, 0, 9, 1) + 5f(0, 13, 14, 3) 


where a, b, c, d, e, f are any nonnegative numbers. 
Setting a = ky, b = ko, 5c = ka, Sd = kẹ, Se = ks, Sf = ke, we have 
for X the representation 


X =k, (5, 0, 3, 0) + k2(0, 5, 4, 0) + k3(3, 0, 2, 0) + 
+ k4(0, 3, 3, 0) + ks(13, 0, 9, 1) + k6(0, 13, 14, 3) (21) 


where ky, ko,., kg are any nonnegative numbers. 

To conclude this point, we note that having proved the existence 
of the fundamental set of solutions for any homogeneous system of 
linear inequalities we have thereby proved Theorem 2’ of Section 7 
on the structure of any convex polyhedral cone. 


11. The Solution of a Nonhomogeneous 
System of Inequalities 


Now that we have learnt how to construct the general solution of 
a homogeneous system of inequalities it won't be difficult to solve 
a similar problem for an arbitrary, i.e. nonhomogeneous, system of 
inequalities. 
Let 
aixi T d3X3 +. + A,X, Z 
X, 2 


bixi + bax2 +. + bx, (1) 


€1X4 + €3X2 b CnXn zc 
be an arbitrary system of linear inequalities in n unknowns 
Xis X34 x, Along with it consider the homogeneous system 
33. 
dX, + d2X2 + oe + a,X, — at > 0 
bixi +baxz + = + bx, — bt>0 Q) 
C4X1 + €aX3 +... 6x, — ct 2.0 
in n+1 unknowns Xj. X2, X, f. ’ 
There is a relation of certain kind between solutions of the systems 
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(1) and (2) in that if xj, x2., X„ t is a solution of the system 
(2), with t> 0, then the numbers 


M 


RE deer (3) 
will constitute a solution of the system (1). 
Indeed, consider, for example, the first inequality of (2). It holds 
for the numbers x4, x;,.., Kab Le, 
X; + d5X5 d... + a,X, > at 
On dividing both sides of the inequality by a positive t we have 
4X, + ay +. +a, >a 
but this means that the set of the numbers X,, Xj... X, satis- 
fies the first inequality of (1). A similar argument goes through 
for any inequality of (2), whence our statement. 

It is easily seen that any solution of the System (1) can be 
obtained by the above method. Indeed, if 3, Xa, X, is some 
solution of the System (1) then the numbers Xem. Xa mE 
X, = X, t = 1 satisfy the system (2) and at the same time equalities (3) 
are valid. 

Thus, in order to find all the solutions X,, X,..., x, of the 
system. (1) one should find all the solutions of the system (2) 
o Which £> 0, and transform each of them according to formulas 

Example. Suppose it is required to find all the solutions of the 
following system 

— 3x1 — 4x, + 5x3 > 6 
2x, + 3x2 — 3x3 > —1 


0 (4) 
0 
0 


M 
VV wv 


x3 


Proceeding as pointed out above, write down the following 
auxiliary homogeneous system 


735353 4x; + 5x54 60 
2x, + 3x2 — 3x5 + t>0 
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Adjoin to this system the inequality t > (for we are interested 
only in such solutions for which t 0) to get the system (19), (20) 
of Section 10, the only difference being that we write t instead of 
X4. The solution set of this system is, as shown in Section 10 (see 
(21) of Section 10), given by the formulas 


xı = Sk, +3k3 + 13k; 

x2 = 5k; + 3kq + 13ko 

x3 = 3k, + Ak 4-2ks + 3k, +95 + 14k 
t=ks + 3k, 


Where k,, kj, kg are any nonnegative numbers. Since we are 
interested only in solutions for which t > 0, it should be assumed that 
at least one of the numbers ks, kg is nonzero (strictly positive). We now 
find the general solution of the system (4) from the formulas 


Lo Sk, + 3ks + 13ks 
xi = 3 ——À. 


ks + 3kg 
_ Ska + 3k, + 13k, 5 
a er 6) 
x: 3k, + 4k; 2k, 3k, + 9ks + 14ko _ 
371 ks + 3ks 


It should be emphasized once again that in these formulas K;, 
ka., kg are any nonnegative numbers, at least one of the numbers 
ks, kg being other than zero. 

Having established that solving the system (1) by the above 
method reduces to solving the homogeneous system (2) we have prov- 
€d Theorem 3 of Section 7 on the structure of any convex polyhedral 
region. This may be illustrated by the example of the system (4). Set 


ee e EDEN 


since the numbers kı, kz, ka, ką are arbitrary nonnegative, so are 
the numbers ki, k, k, k. Further set 
ks 
pis wk 
ks ks + 3kg 


are nonnegative and constrained by the 
ow be written down in the 


the numbers ké and ké 
Condition kj + kg = 1. Formulas (5) can n 
form of a single equality 


Qu. X» X) = KL, 0, 3) +430, 5, 4) + 3B, 0, 2) + 
Ho 13. M 6 
KL, 3, 3) + kz(13, 0, 9) + (0. 3 5) (6) 
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We introduce the notation: 
X,=(5. 0, 3), X2=(, 5,4), X,—-(3,0 2) 


X4= 033, ¥5=(1.09, x, e (o D, 14) 
Bearing in mind the above restrictions on kí, kz, kg, ká as well 
as on kg, kg, equality (6) can now be interpreted as follows: 
the solution set of the system (4) is (X1, X5, X, X4) + (X5, Xs». 
Thereby the statement of Theorem 3 of Section 7 is proved for 
the system (4). 


12. A Linear Programming Problem 


Linear programming is a relatively new field of applied mathema- 
tics which developed in the 1940s or 1950s in connection with 
the tackling of varied economic problems. 

As a rule problems occurring in economics, and in production 
planning in particular, are problems of finding the most profitable 
variant. Each of them calls for an answer to the question as to how to 
use the limited available resources to obtain the utmost effect. Until 
recently the only method of solving such problems has been by ordi- 
nary rough calculation, estimating "by sight", or else by looking 
all the Possible variants over to find the best. The situation is 
changing now, Over the past few decades the complexity of produc- 


xamples of linear programming problems. 
turning out items of two types has an 


the first type and so many items of the Second per diem) which 
would ensure the largest profit for the factory, 
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The desired production plan is specified by means of two nonne- 
gative integers x, y (x being the number of items of the first type and 
y that of the second) which must satisfy the following conditions *: 


- 3x + y < 300; x+y<150 
2x+y is maximal 


In other words, from among the nonnegative integer solutions 
of the system 
3x + y < 300 
s (4) 
x+y< 150 
one should pick the one imparting the largest value to the linear 
function 


f=2x+y 


Given an xOy rectangular coordinate system, the solution set of the 
System (1) will be represented by the shaded polygon of Fig. 47. 


* ondition comes from the assembly shop. Indeed, it can turn 
out ue hore of the second type instead of one of the first. Hence, in 
terms of second-type items the shop's entire production is 3x 4- y articles; 


this number must not exceed 300. 
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From the same figure one can see that the point P (75, 75), one 
of the vertices of the polygon, is the solution of the problem. . 
Indeed, consider a straight line 2x+y=c where € is a certain 
“number. Denote the line by l.. As the number c increases, the line 
l, shifts “upwards” (at the same time remaining parallel to its initial 
position) The largest value of c at which the line /, still has 
points in common with the shaded polygon is that value ofcat which 
the line passes through the point P: Hence at this point the function 
2x + y attains its largest value (in comparison with its values at the 
rest of the points of the polygon). yer. dM. k 
The example discussed is certainly very primitive, still it gives an idea 
of the nature of linear programming problems. In all of them it is 


required tô find the maximal (or minimal) value of some linear function 
of n variables 


f=cCyx, + 2X2 tence CaXR 


on condition that all these variables obey a given system of li- 
near inequalities (which must without fail contain the nonnegativity 
constraints on the variables: x; > 0...., x, 2 0). In some problems (such 
as in the example above) an additional requirement is imposed on the 
variables, that they should be integers. 

Example 2 (diet problem). There are several kinds of food at one’s 
disposal. It is necessary to compose a diet which would, on the one 
hand, satisfy the minimal human body requirements for nutrients 


(proteins, fats, carbohydrates, vitamins, etc.) and involve the smallest 
expenditures, on the other. 


Consider a simple mathematical model of this problem. 
Suppose there are two kinds of food, F, and F;, which contain 
the nutrients A, B, C. It is known how much nutrient of one or another 


kind is contained in 1 Ib of F, or F;; this information is presented 
in the following table. 


A B C 
F,11b ay b, cy 
F, 1 Ib az b; C2 


In addition to these data we kno 
of human body for A, B, C (res 
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W that a, b, c are daily requirements 
pectively) and that s, s; are costs 


of 1 Ib of F, and F; (respectively). It is necessary to compute the 
amount x, of the food F, and the amount x; of the food F so that 
the required amounts of the nutrients should be ensured with minimal 
expenditures on the food. 

Obviously the total cost of the food will be 


S = sx; d 52X2 


The total amount of the nutrient A in both kinds of food is 
1X4 + aX. It must not be less than a: 


dX; + ü2X3 > A 


Similar inequalities must hold for B and C: b,x, + b2x22b, 
CiX1 +x; >c. Thus we arrive at the following problem. 
We are given the system 


aX; d3X532a 
bx, + box. > b (2) 
CX, T C2X22 C 


of three linear inequalities in two unknowns X;, X2 and the linear 
function 


S=s 1X; +S2X2 


It is required to pick from among the nonnegative solutions (x;, 
X3) of system (2) such a solution that the function S attains the 
smallest value'(is minimized). Various problems can be reduced to 
such schemes, viz. alloy problems and problems in fuel-oil blending, 
feed mix problems and problems in fertilizer mixing, etc. 
_ Example 3 (transportation problem). Coal mined in several deposits 
is shipped to a number of consumers, to factories, power stations, 
etc. It is known how much coal is mined at each of the deposits, 
say, per month and how much is required by any of the consumers 
for the same term. One knows the distances between the deposits and 
the consumers as well as the traffic conditions between them; consi- 
dering these data one can calculate the cost of transporting each ton 
of coal from any deposit to any consumer. It is required that the 
transportation of coal should be planned in such a way under these 
conditions that the expenditures on it shall be minimal. 
Assume for simplicity that there are only two deposits Dı, D; and 
three consumers Ci, Ca, C3. The amounts of coal at D, and 
-D, are a, and a; respectively; let the requirements of Cy, C2, C3 
be bı, b», bs respectively. Let us suppose that the total stocks of 


87 


coal are equal to the total requirements: 
a, + a5 =b, +b, +b, 
such an assumption is quite natural. Finally, the numbers c; (i — 1, 


2; j — 1, 2, 3) are given, denoting the costs of transporting a ton of 
coal from D, to C; The problem is to find the six numbers 


Xin X12, X13, X21, X22, X23 


where x, is the amount of coal due to be shipped from D; to C, 
For the convenience of reviewing we make the following table: 


Amount 
shipped 


From D, x 3 ay 


From D, 2 x KEJ 


Amount 
delivered 


The total amount of coal shipped from D, must be 4,; hence 
we have the condition 


X11 tX12+X13 =a, 
A similar condition must hold for D,: 
X21 1:X22.-- X38 5-05 
The total amount of coal delivered to C, must be b,; hence 
X11 c X4 =b; 
Similarly we Bet the conditions 
Xiz-X;3-b;, x13 X33 =b; 


We assume that the cost of transportation is dire 
to the amount of coal conveyed, i.e. transportati 


D, to C, costs cx; The total cost of transportat 


S = 6X1 + C1212 + 13X13 + cy, 


ctly proportional 
on of coal from 
ion will then be 


X21 + €22X22 + C33X23 


Thus we arrive at the following problem. 
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Given the system 
X11 + X45 X13 = 44 
X21 + X22 + X23 = d5 
Xi + X21 =b; (3) 
X12 + X22 = b; 
X13 + X23 = b3 


of five linear equations in six unknowns and the linear function 
S, it is required to pick from among the nonnegative solutions 


X11, X12, X13, X21, X22; X23 


of the system (3) such a solution that the function S attains its 
smallest value (is minimized). 

The problem can certainly be formulated in a more general 
form, i.e. with any number of deposits and consumers. It was given 
the name of “transportation problem" and was one of the first to be 
Successfully solved using linear programming methods. 

We have considered in a simplified form three linear programming 
problems. Despite the diversity of their subjects the problems have 
much in common in their formulations. In each problem one seeks 
the values of several unknowns, it being required: 

(i) that these values should be nonnegative; 

(ii) that these values should satisfy a certain system of linear equa- 
lions or linear inequalities; 

(ii) that at these values a certain linear function should attain 
à minimum (or maximum). 

We remark that condition (ii) seems to be a disjoining rather than 
à uniting condition, for in one case the unknowns must satisfy equa- 
tions and in the other they must satisfy inequalities. But we shall sce 
later that one case is easily reduced to the other. 

It is problems of this kind that linear programming is concerned 
with. More strictly, linear programming is à branch of mathematics 
which studies methods of finding the minimal (or maximal) value of a 
linear function of several variables provided these satisfy a finite num- 
ber of linear equations and inequalities. The number of variables 
and that of conditions (equations or inequalities) may of course 
be arbitrary. In actual problems these numbers may be very large 
(some ten or several tens or even more). 

Let us put the above verbal statement down in strict terms of 
formulas. The general mathematical statement of a linear prog- 
ramming problem looks as follows. 
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Given a system of linear equations 
41X1 + d15X5 E dix, = b, 
d51X1 + d25X3 + ... + d55X, = by 


cide OW ‘GLA TOE) dc Rae ERE: (D) 
Um1X1 + Am2X2 +... + AmnXn = Dy, 
and a linear function 
J= CX, + C22 Lt ex, (1) 
Find a nonnegative solution 
EXE Oh o Okan ue 0 (II) 


of the system (I) such that the function f assumes the minimal value (is 
minimized). 

The equations of (I) are called the constraints of a given problem. 
Strictly speaking, conditions (III) are also constraints; it is not 
customary, however, to call them so since they are not characteris- 
tic of a given problem alone, but are common to all linear program- 
ming problems. 

Of the three examples considered above only the—last (the 
transportation problem) corresponds, one would think, to the state- 
ment just given. The remaining two look slightly different since 
all the constraints they contain have the form of inequalities and 
not of equations. It is possible, however, using a simple method. 


to make inequality constraints 80 over into the equivalent constraints 
given in the form of equations. 


Assume, indeed, that the constraints of a given problem contain 
the inequality 


4X; FaXjct.ctax-bm-0 (4) 


We introduce a new, so-called additional, unknown X44; Connect- 


ed with the unknowns Xi X2,..., X, by the equation 
Xi + d3X3 Fax, +b =x 


n+l 

Inequality (4) is obviously equivalent to the nonnegativity condi- 
tion on x,,,. If an additional unknown is introduced for each 
of the inequalities contained in the system of constraints for a given 
problem, requiring in addition that all additional unknowns should 
be nonnegative, the problem will assume standard form (I), (II). 
(III) although involving a larger number of unknowns. We demon- 
strate this method using the diet problem as an example. The 
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system of constraints consists of three inequalities in this problem, 
so we introduce three additional unknowns xs, X4, Xs. 
As a result the constraints assume the form of equations: 


01X4 .d2X2 — X3 = 4 
bixi + b2x2 — x4 — b 
C1X, c C2X2 — X5 =C 


It is required to find among all the nonnegative solutions of this 
System such a solution which minimizes the function 


S = S1X1 + 52X2 


Certainly, what we are essentially concerned with in the required 
solution are the values of x, and x3. 

It is frequent for a linear programming problem to require that 
the maximum, and not the minimum, of a linear function f should 
be found out. Since max f= — min (— f) (establish this for yourself), 
one problem is reduced to the other if f is replaced by — f. 

We have shown that a linear programming problem having inequa- 
lity constraints can be reduced to a problem with equation constraints. 
It is highly interesting that the converse should also be possible, i.e. 
any linear programming problem can be stated in such a way that 
all constraints will assume the form of inequalities. We shall not 
dwell on this, however. 

In conclusion we shall make a few remarks concerning the adopted 
terminology. Any nonnegative solution of a system of constraints 
is called feasible. A feasible solution yielding the minimum of a 
function f is called optimal. It is the finding of an optimal solution 
that is our goal. The optimal solution, if any, is not necessarily unique; 
cases are possible where there may be an infinite number of optimal 
solutions. 


13. The Simplex Method 


Actual linear programming problems tend to contain a large 
number of constraints and unknowns. It is natural that solving such 
problems involves a large number of calculations. This difficulty is 
overcome with the help of high-speed computers. The algorithm a 
computer programme is based on may be connected with a specific 
class of problems. Thus, for example, there are very simple algorithms 
for solving the transportation problem which are conditioned by 
the peculiarities of its system of constraints. There exist, however, 
general methods as well which allow a solution of a linear prog- 
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ramming problem to be found in a visible number of steps. The 
first and foremost of them is the so-called simplex method with some 
of its modifications *. ; 

1°. A description of the simplex method. So let us consider a 
linear programming problem. Given a certain system of linear equa- 
tions in n unknowns x), x5... X, and a certain linear function Í it 
is required to find among the nonnegative solutions of the given 
system such a solution which minimizes the function f. 

To begin working by the simplex method it is necessary that thc 
given system of equations should be reduced to a form such that 
some r unknowns are expressed in terms of the rest, the absolute 
terms in the expressions being nonnegative. Assume, for example. 
that n=5 and that the unknowns expressed in terms of the rest 
are Xi, X2, x3. Hence the system of constraints is reduced to the 
form 

X, = X + 4X4 + OsXs 
x2 = B + Bax, + Bsxs (1) 
X3 = Y + Y4X4 + Y5X5 
where 
a>0, B>0, y>0 (2) 


Is it really possible to reduce the system to such a form and how 
is this to be done? This question will be considered later on (see 
sect. 2°). The unknowns Xis X2, X3 in the left members of the 
system (1) are called basic and the entire set (Xj, X5, x3} which 
for brevity we shall denote by a single letter B is called a 
basis; the remaining unknowns are called nonbasic or free**, Sub- 
stituting into the original expression for the function f the expres- 
sions in terms of the nonbasic unknowns of (1) for the basic un- 
knowns, we can write down the function f itself in the nonbasic 
unknowns x4, x5: 


fet eaxa + esxs 
We set the nonbasic unknowns equal to zero 


X4=0, x,-20 


* The name of the simplex method has nothing to do with the essence 
of the matter and is due to a casual circumstance. 
** The latter name is due to the fact that when finding the solutions of 


the system (1) (irrespective of a linear programming problem) these unknowns, 
may be assigned any values. 
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and find from the system (1) the values of the basic unknowns: 
4 =O, x2—= 8, =y 
The resulting solution of the system 
(a, B. y, 0, 0) 


will, due to (2), be feasible. It is called a basic solution corresponding 
to the basis B = (xi, X2, X3}. For the basic solution the value of the 
function f is 
UL 

The solution of a linear programming problem by the simplex 
method falls into a number of steps. Each step consists in going over 
from a given basis B to another basis B' so that the value of f should 
decrease or at least not increase: fy € fj. The new basis B’ is derived 
from the old basis B in a very simple way: one of the unknowns is 
eliminated from B and another (a former nonbasic) unknown intro- 
duced instead. A change in the basis involves a corresponding 
change in the structure of the system (1), of course. After a certain 
number k of such steps we either come to a basis B for 
which fy) is the desired minimum of the function f and the cor- 
responding basic solution is optimal, or else find out that the 
problem has no solution. 

We illustrate the simplex method by some examples. 

Example |. Let a given system of constraints and a function f 
be reduced to the following form: 


Xy=1— x,-42x5 
x =2 +2x4— Xs 


X423 —3x4 —. Xs 


jJ — X4 — Xs 


Here the unknowns x;, X2, x3 form the basis. The corresponding 
basic solution is 


(1, 2, 3, 0, 0) 


the value of f is 0 for this solution. 

We find out whether the solution is optimal. Since xs has in fa 
negative coefficient, we may try to decrease the value of f by 
increasing xs (while retaining a zero value for x4). Care should 
be taken in doing so, however, since changing xs will lead to changes 
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in the values of x,, x2, x3 and it is necessary to see to it that none 
of them become negative. ; 

There is no such danger for x,, since increasing x; results 
in the increase of x,. Examining x, and X3 we find that xs can 
be increased to 2 (at most, otherwise X3 will become negative) which 
yields x, — 5, x; —0, x, — 1. As a result we have a new feasible I 
solution (5, 0, 1, 0, 2) in which the number of positive unknowns 
is equal to three as before. The value of f for this solution is 
equal to — 2. 

The new basis now consists of Xj, Xs, X3. To make appropriate 
modifications in the system of constraints, it is necessary to express 
these unknowns in terms of x, X4. We begin with the equation 
for x2 (the new nonbasic unknown) solving it for xs (the new basic 
unknown): | 


Xs = 2 + 2x4 — x5 
we then express x,, x3 and f: 


SUSU x4 -+2(2 2x4 — xj) 


U ES X4— (2+ 2x4 — xj) 
Thus the problem is reduced to the form 
X1 = 5 + 3x4 — 2x, 
X5=2 +2x4— x, 
x3 = 1 — 5x4 + X3 
J= 2 Kat Xp 
The new basic solution is 
(5, 0, 1, 0, 2) 


the value of f for this solution is equal to — 2. This completes the 
first step of the process. 

Let us see if it is possible to decrease the value of f still 
further. The coefficient of X4 in the expression for f is negative, 
therefore it is possible to try to decrease f by increasing x, (without 
changing x; — 0) There is nothing in the first and the second 
equation to prevent it. and it is seen from the third that x4 can 
be increased to 1/5 (at most, for otherwise X3 will become negative). 
Setting x, = 1/5, x; — 0, we have x, = 28/5, x, = 12/5, x4—0. Asa 
result we arrive at a new feasible solution (28/5, 0, 0, 1/5, 12/5). 

The new basis now consists of x,, x.. x4. The equation for x3 (the 
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new nonbasic unknown) is solved for x4 (the new basic unknown): 
1 1 1 
X4=5 3*3 EE 3*2 
and the resulting expression is substituted into the other equations. 
As a result the system assumes the form 


uid 


The new basic solution is 


(28/5, 0, 0, 1/5, 12/5) 


the corresponding value of f being equal to — 11/5. This completes 
the second step of the process. 

Since in the last expression for the function f both unknowns 
xs and x, have positive coefficients, the minimum of f is attained 
when x4 = x; = 0. This means that the last basic solution (28/5, 0, 
0, 1/5, 12/5) is optimal and the required minimum f is equal to 
— 11/5; thus the problem is solved. 

In the example analysed the process has terminated in finding 
an optimal solution. Yet another termination of the process is possible. 
To illustrate it, we solve the following example. 


Example 2. 
xy = 4X 
x4 =2— xı + 2X2 
f2*,—X2 


Here the basis is formed by the unknowns xs, xa. The basic 
solution is of the form 


. (0, 0, 1, 2) 


the corresponding value of the function f being equal to 0. 
In the expression for f the coefficient of x, is negative, therefore 
we try to increase x; (without changing x; — 0). The first equation 
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does not prevent us from doing so, and it is seen from the second 
equation that x, can be increased only to 2, which yields x3 =3, 
X4 —0. The new feasible solution is (2, 0, 3, 0); for this solution 
the value of f is equal to — 2. 

The new basis now consists of x3, Xi. We solve the second 
equation for x, and substitute the resulting expression into the first 
equation and f. The problem assumes the following form 


X= JioÀx cx 
= 2—Xx4-2x; 
J= —2 + x, — 3x: 


The new basic solution is 
(2, 0, 3, 0) 


the Corresponding value of f being equal to — 2. 

In the last expression for f the coefficient of X2 is negative. We 
want to see how much we can decrease the value of f by increasing 
X2 (while keeping x4 = 0). For this Purpose we look over the x; terms 
in both equations of the System and notice that both coefficients 
of x; are positive. Thus X? can be increased without limit while 
keeping x; and X1 positive, f taking arbitrarily large absolute nega- 
tive values. So min J=- co (it is said that the function f is not 
bounded from below) and there exists no optimal solution. 

As is seen from the above examples, each successive step in the 
simplex method is based on the choice of a negative coefficient for 
some nonbasic unknown X, in the expression for f. If it turns out 
that there are several negative coefficients, one m 


X; exchange roles with X, i.e. to make 3 A i : 2 


| 


simplex method. We required as a preliminary condition that the Sys- 
tem of constraints should be reduced to the form * 


Xp SOHO Xp Ho FOX, | 
x2 = p +t Biss Mra the F Boss \ (3) 


= eagle oe IRE TEE UTR SS eos ls f 


X EY YeaXac F Yota | 


where à > 0, B > 0, . y 20: then we say that the unknowns x, 
X2, ... X, form the basis. 

In many linear programming problems the basis can be directly 
perceived, and in others it has to be found. We shall consider one 
of the methods of finding the basis generally known -as “the method 
of an artificial basis”. 

We shall analyse an example to illustrate the method. 

Example. Given the system of constraints 


2x, — 2x; + 2x3 — xy — Tx; = — 
=x; + 3X2 — 2x3 + x4 + 3x5 =S 
it is required to solve it for a certain original basis. 
Let us first of all transform the system so that the absolute 
terms of the equations should be nonnegative. To do this we should 


multiply both members of the first equation by — 1. We 
get the following system 


O — 2xs + Xe + Ts = s (4) 
—x, + 3x2 — 2x3 + X4 3x4 =5 


We now introduce auxiliary, or artificial, unknowns yı, y». (one 
for each equality) in the following way: 


yı = 4 — (2x, + 2x2 — 2x3 t x4 + m 


yo = 5—(— xi 3x2 — 2X3 + x4 + 3x5) 6 


, Obviously all solutions of the original system (4) can be obtained 
by taking all the solutions of the system (5) which satisfy the 
conditions y; — 0, y; — 0, and keeping only the values of x4... x5 
in each of those solutions. 

The unknowns y, y; form the basis in the systém (5). Suppose 
that starting from this basis we succeed in going over to another 

* Or to a similar form where the left members contain some other r 
unknowns rather than xy, Xz x,; the remaining n — r unknowns must 
enter into the right members. 
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basis which contains no artificial unknowns. Then omitting the terms 
containing y, and y; from the resulting equations (or, equivalently. 
setting y, =0, y; — 0) we get a system equivalent to the original 
System (4), this system being solved for a certain basis. 

It remains to decide how to transform the unknowns Vis ya of 
the system (5) into nonbasic unknowns. It is remarkable that we 
should be able to use the simplex method for this purpose too; 
namely, using this method we shall work out the problem of mini- 
mizing the function 


F= y, +y 


under constraints (5)and the conditions X120,.,x5 20, 9,20, y; 20. 
Here the premises under which the simplex method begins to work 
are satisfied since the system (5) is of the required form (it has 
the original basis Yi, Y2). In some steps we arrive at the desired 
minimum. Since F > 0 and hence min F Z 0, two cases are possible: 

l. min F 7 0. This means that the System (5) has no nonnegative 
solutions for which y, — 0, Y2 — 0 (otherwise min F'— 0). Therefore 
the original system (4) will have no nonnegative solutions in this case. 
Hence it will follow that any linear programming problem with such — 
a system of constraints is insoluble. 

2. min F = 0. Suppose that in this case (xf, ..., x8, y9, y9) is an 
optimal solution. Since y? + y2 = min F = 0, we have y = 0, y2 =0. 
It follows that (x9..... x5) will be a nonnegative solution of the original 
system (4). 

Thus in the case of min F — 0 the System of constraints (4) will have 
at least one nonnegative solution. | 

Further, if on accomplishing the process (i.e, when min F — 0 has 
been attained) it is found that the artificial unknowns Vis yi are 
among the nonbasic ones, we have achieved our object —we have | 
separated the basis from the unknowns X1. Xs. If however, it becomes 
clear that some of the unknowns Yı» ¥2 Still belong to the basis, then | 
some further transformations are required, 

So we solve the problem of minimizing the function 


Foy +y = I4 (xx 2x3 + x4 + 7x5)] + | 
*[S5—-(— x So 2x3 + x4 + 3x5)] = 
=9 + 3x; — 5x. + 4x, — 2x4 — 10x, 


In the expression for F the coefficient of x, is negative, we there- — | 
fore try to decrease the value of F by increasing X4. Comparing the 
Coefficients of x, and the absolute terms in the equations of (5) wen 
find that x, can be increased only to 4, y, then vanishing. The — | 
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new basis consists of x4, y;. We solve the first equation for x, and 

substitute the resulting expression into the second equation and 
into F. The problem takes the form 

X4 = 4 + 2x; — 2x5 2x3 EAEAN 

Yo =1— x,— x + xs + yif 

F21— xj— x; +4x5+2y, 

This completes the first step of the process. As a result we have 
succeeded in making the unknown y, nonbasic. 

In the new expression for F the coefficient of x, is negative. Looking 

at equations (6) we see that x; can be increased only to 1, y2 then 


vanishing. The new basis consists of x4, xı. On transforming equa- 
tions (6) and the expression for F (transforming the latter is now 


Optimal, however) we have 
xa = 6 — 4x2 + 2x3 + Xs + yı — 2y2 
xı =1— x +4xs +y — Y (7) 
Fi Ve 
Thus as a result of the second step the artificial unknowns y1, Y2 have 
become nonbasic. Eliminating from equations (7) the terms having 
Yi, y», we get 


(6) 


x4 = 6 — 4x2 '+ 2x3 + a 
x, =1= X; + 4x5 


This means that we have separated a basis in the original system 
(4). Thus the problem is solved. l í 

As already noted, when applying the simplex method to poine 
the minimum of the function F equal to the sum of the artificia! 


unknowns, it may happen that by the tinens process m omes 
] i A ttained) some o 
When min F..0 haud Be base These unknowns can also 


unknowns still make part o I i 
be transformed into nonbasic ones by simple methods which we 


shall demonstrate using an example. ; 
Suppose, for example, that after a certain number of steps the 


problem has assumed the following form 
3 3 1 


1 
Xi 5X3 — 5X4 — 9*5 73h 


(8) 
Yo = Xi— X3 + ldix,44 5x5 + 2y1 
ys = 6x, + 6x3 + 6x4 4+ 4x5 F2yi - 


F = Tx, + 5x3 + 17x4 + 9x5 + 5y1 
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We see that the function F has already attained the minimum 
equal to zero, it may therefore be omitted from further considera- 
tions. The unknowns y; and y; are still part of the basis. But 
the absolute terms in the equations for y, and y, are zero and not 
by chance, for min F — 0 implies that in the resulting basic solution 
y2=0 and y; — 0. We now make use of this circumstance to make 
y, and y, leave the basis. To do this we observe that there is a 
negative number among the coefficients of the unknowns Xpes X5 
in the equation for y;: it is the coefficient — 1 of X4. We make 
the replacement x5 € y, (we make the unknown X3 basic and make 
y2 become nonbasic). To do this we solve the second of the equations 
for x3 and substitute the resulting expression for x3 into the other 
equations. Since the absolute term in the second equation is zero, as à 
result of such an operation the absolute terms in the equations 
will remain unchanged (and hence nonnegative). 

We get: 


1 5 3 
x2 =— +X, + 15x4 +7x5 + —y, — >y, 
2 2 2 


X= My + Ixy + $x, 2y - y, (9) 

V3 = 12x, + 72x4 + 34x5 + 14y, — 6y5 
Now we have no negative coefficients (of x,..., xs) in the equation 
for ys, therefore we shall not be able to make y, nonbasic. But 
we should not be much annoyed at this. For if some numbers x. 


^» Xs Satisfy the original system of constraints, then y,, y; and ys 
must be zero. Thus ite 


12x, + 72x4 + 34x; = 0 


Since here all the coefficients of the unknowns have the same 
sign, this must yield 
Xi X42xX,—0 


these equalities are thus consequences of the original system of 
constraints and the conditions x, > 0 (i = 1,.., 5). If we regard them 
as satisfied the system (9) reduces to 


D 
N| = 


x3=0 
So subject to the nonnegativity of the unknowns the given system 
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of constraints allows only one solution: 


xu x3=0, x;20, x4=0, x5=0 


14. The Duality Theorem 
in Linear Programming 


There are so-called duality theorems in various branches of mathe- 
matics. Each of them allows the constructing for any statement 
in a given theory, according to a certain standard rule, of another 
statement in such a way that the validity of the first statement 
automatically leads to the validity of the second. A remarkable 
example of a duality theorem is met in linear programming too. 

Suppose a certain linear programming problem is given; we shall 
call it the original problem, or problem A, below. The constraints 
are given in the form of inequalities. So we have the system 


Ay 4X1 + 12X2 + + AinXn + b, >0 
a21X1 + d23X2 + a + aan + b2 20 (1) 


Amt Xi + Am2X2 Foe Xa + by, > 9 
of m linear inequalities in n unknowns and a linear function 
fex + C2X2 c CaXn 


We must pick from all nonnegative solutions (x; > 0, x2 2 0... 
X, > 0) such that the function fattains the largest value (fis maximized). 

We connect with problem A another problem which we call a dual 
problem relative to problem A or problem A’. It is stated as follows: 
given the system 


Ay V1 H danY2 o H gay H Cn «0 
alities in m unknowns and the linear function 


Q = biyi + Days + bus 


of n linear inequ: 


It is necessary to pick all nonnegative solutions of the system (1’) 
such that the function @ attains the smallest value possible («p 
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is minimized). Comparing problems A and A’ we notice: : 

1. That the coefficient of the jth unknown’ in the ith equation 
of the system (1') is the same as that of the ith unknowns in the jth 
equation of the system (1); 

2. That the absolute terms of the inequalities in each of the prob- 
lems coincide with the coefficients of the unknowns in the linear 
function of the other problem; 7 

3. That in the system of inequalities of problem A the inequalities 
are all of the 2 0 type, it being required in the problem that the 
maximum of f should be attained, On the contrary, in the system of 
inequalities of problem A' the inequalities are all of the <0 type, 
but in return it is required in it that the minimum of @ should 
be attained. 

One of the main theorems in linear programming, the so-called 
duality theorem, states the following. 

Duality theorem. If an original problem is solvable, so is its dual, 
the maximum of the function f being equal to the minimum of the 
function ọ: 


max f=ming 


We shall prove this theorem by reducing it to the question 
of compatibility of a certain System of inequalities. 


A To make the proof more convenient to follow, we break it down 
into several stages. 


Stage 1. Lemma. If x». x9 is some nonnegative solution of the 
System (1) and y... Y» is some nonnegative solution of the system 


(1^, then for these solutions the values of the functions f and Q are 
connected by the inequality 


Jo & Qo 


Proof. We consider the inequalities of the System (1) where the 
values x?...., x? are substituted for Xi X, We multiply the first of 


the inequalities by y?. the second by y$, etc., and then add all the 
Inequalities obtained: 


(a11 yix + + auyox?) + by y? pou b,ys > O 


(one should bear in mind that we are multiplying the inequalities 
by nonnegative numbers, therefore the signs of the inequalities 
remain unchanged). In the same way, we multiply the first inequal- 
ity of the system (1’) by x9, the second by x9, etc., and then add 
the resulting inequalities together: 


00 0 ^ 
(ay x1 +.. + a, VX) + CX) +... $6,x°< 0 


* 


In both cases the brackets contain an expression equal to a sum ` 
of terms of the form ajy?x$ over all i= 1.., m, j= 1.., n. Hence 
the two expressions in brackets coincide. But then 


c x9 +. exo < Diy? 4-4 b,yo 
Or f; € Po. Thus the lemma is proved. 
Stage 2. Reducing problems A and A' to the solution of a certain 


System of inequalities. 
Consider the following “combined” system. of inequalities: 


d11X1 T + Xn 


F (S) 
. iyi + oe pays + Cn SO 
CyXy F os F Can | —biyi — -~ — bmYm 20 

It is seen that it is made up of the system (1), the system 


(l') and the inequality f— 20. The unknowns of the system 
(S) ‘are Soe ere Wivam- du^ (there are n+m unknowns in all). 


We first of all establish the following fact. 

If the system (S) has nonnegative solution X1... x9, 0... Voa then 
the number: xQ,., x9 give a solution to problem A and the numbers 
yt, y° a solution to problem A’, with fo = Go. 1 

We'shall linger on this proposition a little longer to emphasize 
the principal role it plays. What is remarkable about it 1s that 
a linear programming problem, i.e. a maximization problem, reduces to 
Solving a certain system of linear inequalities without any maximiza- 
tion requirements. In fact, the solution of the system (S) (in the 
Tange of nonnegative values of unknowns) is of course not a bit 
easier than the solution of the original linear programming problem 
(problem A); however, the very possibility of such a reduction is 
Very curious. 


Now we prove the above statement. 
the numbers x?...., x? are nonnegative an 
larly the numbers y}, y? are nonnegativ 
Over, these numbers satisfy the inequality 


Jo > Qo 


First of all, it is clear that 
d satisfy the system (1); simi- 
e and satisfy (1). More- 
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(following from the last inequality of the System (S)) On the 
other hand, by lemma we have 


Jo & Yo 


Hence, f; = Qo. 1 ; . 
Further, if x..., x, is some nonnegative solution of the system (1), 
then by the lemma we again have 


f< 9o 


Comparing this and fo =Po we get f< fo, whence it follows that 
fo is the maximal value of f. 


Similarly, if yj... y. is some nonnegative solution of the system 
(1), then by the lemma we have 


foo 


comparing this and f, = Po We get Qo <Q, i.e. Po is the minimal 
value of q. This proves the above proposition. 

Stage 3. Completing the proof. Now it remains for us to show 
the following: if problem 4 has a solution, then the system (S) has a 
nonnegative solution. For then, as shown above, fo =o, ie. 
max f— min ọ 


of the <0 type 


! - But this is easy to mend by writing the system 
(S) in the form 


1X4 + i 1X, +b,>0 
[LET tb.20 
E (s) 
7ny ys — C1, > 0 
= Gini — + — AmnYin — Cy > 0 
C€pXi-b s rex, |S by, —.— bn 20 


So suppose that the system (S') has no nonnegative solutions. 
According to the corollary of the theorem on incompatible systems, 
there will occur nonnegative numbers Kis s his los (mnl 
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numbers all in all) such that* 
aiki + + auk + C18 S 9| 


à da sn erae M uie d o ac Í (2) 
An ky + oe + AnKa + CnS < OF 

— a;l inla — bys <0 / 

sce 10 LEE E (2) 
Anil "ey Ganla = bys < 0 

bi ky Fio bue, elt. = Culp <0 (3) 


We first of all show that the number s is nonzero. Indeed, suppose 
this is not the case, ie. that s=0. Consider some nonnegative 
solution x$..., x° of the system (1). Arguing in the same way as we 
did in proving the lemma, we find that 


(ay Kx +... + OyankinX®) + biki + + bak > 0 
But the parenthesis < 0 (this results if we multiply the first inequality 


of (2) by x9, the second by x? and so on, bearing in mind s — 0, and 
then carry out addition). Hence 


biki +... + bmKm > 0 (4) 
after which it follows from (3) that 

eil + -e + Cal > 0 
Further, by multiplying inequalities (2’) (where s = 0) by an arbitrary 


Positive number à we get 
aidli +. + aA, > 0 


amidli + agg, > O 
If we now add to each inequality of (1), where Xis X, are 
substituted for by x0... Xm the corresponding inequality of (5), we 
find that the numbers 
xD Alas x9 + Al, 


(5) 


S—— 

* The numbers kis- Kg» lis lpp s are precisely those numbers by which 

We multiply the first, second, ... (m+n 4- 1)th inequalities of the system 

) to get (after addition) the inconsistent inequality a,X t + [Ru 

T a y, +. + ds ya + d > 0, where aj... d, dj... aj, € 0 and d is a negative 
number, 
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also constitute a nonnegative solution of the system (1). For this 
solution the value of the function f is 


recy fares Cae (Cyl, F eek Caly) 


and since the parenthesis is strictly positive, the value of f increases 
indefinitely with A. This means that max f= æ, i.e. that problem A, 


Thus s is nonzero. It then follows from (2) that the numbers 
ki/s., k,/S are a nonnegative solution of the system (1), from 
(2) that the numbers h/s.. L/s are a nonnegative solution of 
the system (1), and from (3) that for these solutions 9—/f «0. 
But this is contrary to the lemma. So on supposing that the system (S) 


Example. Find the maximal value of the function 
f 2x2 + 12x, 
provided the variables x,, x5, x, are nonnegative and satisfy the 
inequalities 
X(—X;— x4,4220 
—X;,—xX;—4x,4-120 
Solution. Let us call the set problem problem A. The dual problem 


(problem A') must be stated as follows: find the minimal value of 
the function 


922v +y 
provided the variables y,. y; are nonnegative and satisfy the inequali- 
ties 
Ji dy <0 
—Jy1— y2+2<0 (6) 
=y: — 4y2+12 <0 


Problem A’ can be solved graphically by representing in the 
Y10y; coordinate plane the feasible region of the system (6). This 
is done in Fig. 48. It is seen from the same figure that the 
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Fig. 48 


ERN Q attains its smallest value at the point (0,3), one of 
x vertices of the region. This value is equal to 3. By the duality 
theorem the maximum of the function f must also be equal to 3. 


15. Transportation Problem 


In Section 12 we have considered a number of specific linear prog- 
ramming problems. Of special interest among them is the transpor- 
tation problem, above all due to its practical significance. There 
is voluminous literature devoted to this and similar problems. One 
must say that the methods of solving the transportation problem 
arerather instructive. They demonstrate the indubitable fact that when 
Solving any special class of problems general methods (say the 
Simplex method) should be employed with caution, i.e. one should 
consider in full measure the peculiarities of the given class of 
Problems. 


e remind the reader how the transpor- 


1°. Stating the problem. Wi 
he general form. There are some, say 


tation problem is stated in t 

m, source points (suppliers) 

Ai, A25. An 

n, destination points (consumers) 
Bi, By,..., B, ‘ 


We use a; to denote the quantity of goods (say in tons) concentrated 
at point Ais], Ze: PP and b, to denote the quantity of goods 
i 


and so many, say 
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expected at point B; (j= 1, 2.., n). We assume the condition 
d; +a, + + dy =b; b + +b, 


implying that the total stock of goods is equal to the summed 
demand for it. In addition to the numbers a;, b;, we are also given 
quantities c;; denoting the cost of transporting a ton of 
goods from point A; to point B;. It is required to develop 
an optimal shipping schedule, i.e. to compute what quantity of goods 
must be shipped from each source point to each destination point 
for the total cost of shipments to be minimal. 

So our problem has mn nonnegative unknown numbers x;;(i= 1, 
~» m; J = blo n) where x,; is the quantity of goods to be shipped 
from A; to Bj. Lest we should complicate the things with cumber- 
some expressions, we take specific values for m and n, namely, 
m=3 and n=4. So we shall assume that there are three source 
points, A,, 45, Az, and four destination points, Bj, B5, B3, B4. The 
unknown quantities Xi; can be summarized in a table called the 
transportation table: 


Table 1 


Destina- 
tion 
Obr point To By | To B, To By | To B, | Stocks 
point + 
L 
From A, Xin Xia Xi3 Xia ay 
From Az X21 X22 X23 X24. li- a | 
From A; X31 X32 X33 4 X34 | > 88 
Requirements bi bi b, bs 7 
zs 


Since the total weight of the goods shipped from A, to all 
consumption points must be equal to a we have the relation 


Xii + Xiz FX +%44= a; 


Similar relations must hold for the points 4; and 4,. The 
unknowns of our problem must therefore satisfy the equations 


Xii Xi X xad 
X21 X22 FX txa =a 1) 
23 24 2 


Xa, $ X39 X34 +X 


34 — d3 
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But we should.also bear in mind that the total quantity of 
goods delivered to B, from all source points must be equal to by, i.e. 


X11 + X21 + X31 =b; 


Similar relations must hold for the points B5, B3, B4. This leads 
to the equations 


Xp + X21 +%31 =b; 
X12 + X22 + X32 =b2 
X43 + X23 + X33 = b3 
X14 + X24 + X34 = Da 


Notice that equations (1) and (2) are very easy to remember. 
Indeed, the ith equation of the system (1) implies that the sum 
of the unknowns in the ith row of the transportation table is 
equal to a; we shall call equations (1) horizontal for this reason. 
Similarly, the jth equation of the system (2) records the fact that 
the sum of the unknowns in the jth column of the transportation 
table is equal to bj; because of this we shall call the equations 
of the system (2) vertical. 

We ship x, tons of goods from point A, to point B, the cost 
of transporting one ton being cy Hence transportation from A, to 
B, costs cjx, and the total cost of transportation will be 


$2Ecjx, (3) 


(2) 


where the symbol X put before cx; signifies that the quantities 
2, 3, 4 (there 


ty : 3 
CjX, must be summed over all i = 1, 2, 3, and all j = 1, 2 


Will be 12 summands in all). , 
We thus come to the following linear programming problem. 
Given the system of equations (1), (2) and the linear Junction 

(3). Find among the nonnegative solutions of the system such that 

minimizes ction (3). E 
Leelee Mum can be solved by the simplex method 

like any other linear programming problem. Because of the special 

structure of the system of constraints (1), (2) the general procedure 
of the simplex method is, of course, greatly simplified when applied 
to the transportation problem. Here we present a method for solving 
the transportation problem called the method of potentials. It is a vari- 
ant of the simplex method specially adapted for solving the transpor- 


tation problem. — X 
2^. Finding the first basis. As we know, the work by the simplex 
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method is preceded by a preparatory stage, that of finding the first 
basis. In the case of the transportation problem there is a very 
simple and convenient method of finding the first basis called the 
northwest corner rule. The essence of the method is best explained 
by considering a specific example. Suppose we are given three source 
points, 4, A2, As, and four destination points, B,. Bs, B3, By, the 
stocks and requirements being equal to the following quantities: 


a,=60, a,=80, a,=100 
b,=40. b, = 60, b;=80, b,=60 


These data are tabulated in Table 2 below. 


Table 2 
B, B, | B, B, Stocks 
A, 40 60 
4i 80 
4s 100 


Requirements 40 60 80 60 


Let us try to satisfy the requirements of the first destination 
point B, using the stocks of the first source point 4,. This can 
be done in this case since the stocks a, = 60 are greater than the 
requirements b, = 40. 

We therefore enter the number 40 in the cell X11. The require- 
ments of the point b, will be found fully satisfied, and therefore 
the B, column may be temporarily excluded from consideration. It is 
now possible to consider that we have come to new Table 3 
which has three destination points, B}, B}, B4, and three source 
points, A1, Az, Az. the stocks at A, being a, = 60 — 40 = 20, Note 
that in Table 3 the sum of all requirements is as before equal 
to the sum of all stocks. 

We apply the same method to Table 3 and try to satisfy the 
requirements 5; = 60 of the point B, (it assumes the role of the 
first point in Table 3) using the stocks a, = 20 of the point A. It 
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is obvious that we cai i i 
an satish y i i 
SA y these requirements only partially since 


Table 3 
B; B, Bs 
A, 20 20 
15 80 
As 100 
60 80 60 


we opter the number 20, this being the maximum of what can 
of eel gitar from A, to Bp, in the cell x;2. The requirements 
em 2 will be reduced to b; = 40 and the stocks of A; will be found 
i y exhausted. In virtue of this the 4; column of Table 3 may be 
emporarily removed. Thus we come to Table 4 which has only two 
Source points, 4» and 4s. and three destination points, B2, Bz, By. 


now. 


Table 4 
B; Bs Bs 
A; 40 80 
As d 100 
40 80 60 


to reduce successively the resulting 
of all the destination points. 
blem all the stocks of the 


we continue 


In a similar 
tables until we satisfy the requirements 


In virtue of the condition of the prol 
Source points will be found exhausted *. ; 
We get the following values for some of the unknowns in the 


ts, either the row or 


T * When "new" stocks are equal to “new” requiremen 
he column may be excluded at will from the table 
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course of reducing the tables: 


X11 =40, x,,—20, x22=40, x23=40, x33 =40, 
: X34 = 60 (4) 
On entering them in Table 2 we get Table 5. 


Table 5 


Let us agree to refer to those cells in Table 5 which have the 
values of unknowns entered in them as basic, and to the remaining 
cells as vacant. If we consider that the values of the unknowns 
x, which correspond to the vacant cells are zero, then the result- 
ing set of all the unknowns yields the feasible solution of our 
problem. Indeed, it is easy to verify that the sum of the values 
of the unknowns in each row of the table is equal to the 
stocks of the corresponding source point and that in each column, 
to the requirements of the corresponding destination point. Therefore 
equations (1), (2) are satisfied. It now remains to remark that the 
values of all unknowns are nonnegative. 

In the general case, i.e. when there are any number m of source 
points and any number of destination points, the method we have 
described allows us to fill m+n—1 cell in the table. Indeed, in 

“every step we fill exactly one cell, after which one row or one 
column is deleted from the table; there is one exception, where the 
table consists of a single cell on filling which we exclude at once 
both the row and the column. Since the number of all rows is m 
and the number of all columns is n, it is clear that the number 
of all steps and hence the number of filled (basic) cells is m + — 1. 

We now show that the system of constraints (1), (2) can be solved 
for the unknowns corresponding to the basic cells. Hence it will follow 
that these unknowns may be assumed to be basic and the remaining 
unknowns corresponding to the vacant cells in the table to be non- 
basic. : 

To prove this, we join the basic cells with a broken line in the order 
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they arose in the i 
S procedure described above. There will 
broken line of the following form for the example eae 1 


A is not difficult to see now that the unknowns occupying the 
; ee cells can be expressed in terms of the unknowns occupying 
vacant cells. The expressions are found in consecutive order: 


fet for xq, taking the first vertical equation; 

m for x,4, taking the first horizontal equation; 
en for x55, taking the second vertical equation; 
then for x23, taking the second horizontal equation; 
then for x33, taking the third vertical equation; 
finally for x34, taking the third horizontal (or the 


fourth vertical) equation. 


So we have found the basic set for our problem: 


Xib Xi» X22 X25 X33» X34 


This does not complete the solution of the set problem since 
the expressions for the basic unknowns in terms of the nonbasic 
Ones still remain to be found. But these expressions will not be 
required in explicit form at the next stages of solving the transpor- 
tation problem. 

- Solving the problem by the mei 
Of the first basis js but a preparator 
i hen this stage is over, all the unk 

to two groups: 

Xj, basic unknowns, and Xy 
da actual expressions for the basic unknowns in terms of nonbasic 
R es will be needed. As to the function S (the total cost of 
Foon ations it is obligatory to know its expression 1n terms of 

Onbasic unknowns. In other words, it is necessary to express the 
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thod of potentials: The finding 
y stage of solving the problem. 
nowns are found broken down 


nonbasic unknowns 


function S as 
S= a SS (5) 


We shall show how to find coefficients s,, of nonbasic unknowns. 
We assign to each source point A; a certain value (i = 1, 2... m). 
the “potential” of point A, Similarly, to each of the destination 
points B, we assign a value B; (j— 1, 2.., n), the “potential” of 
point Bj. We relate these values as follows: for every basic unknown 

x,, we work out an equation 
9, +B) = Cy (6) 


a 


where c,, is used as before to denote the cost of transporting 
a ton of goods from a point A, to a point B, The set of the equa- 
tions of the form (6) worked out for all basic unknowns x, forms 
a system of linear equations. This system has m+n — | equations 
(as many as there are basic unknowns) and m + n unknowns o, Bı 
(as many as there are source points and destination points taken 
together). It can be shown that this system is always compatible, 
it being possible to arbitrarily specify the value of one of the un- 
knowns and then unambiguously find the values of the remaining 
unknowns from the system. 

Let us fix some one solution (2,..., o. By, f) of the system of 
equations (6) and then calculate the sum 2, + B, for each nonbasic 


unknown x, We denote this sum by chy! 


ba 
tp + Ba = Cp, 
and call it indirect cost (in contrast to the real cost C). It then 


. turns out that in the expression (5) that interests us the coefficients 
of nonbasic unknowns are equal to 


BL b: 
Spa = pq — Cpg (7) 


Formula (7) will be proved in the next section. 

If all the variables Sp are nonnegative, then the initial basic 
solution will be optimal. If, however, there are negative variables, 
Say Spgs among them, then we go over to the next basis. It will 
be recalled that this step begins with reasoning as follows. We shall 
increase X,, (keeping the other nonbasic unknowns zero), If in 
doing so we reach a moment when one of the basic unknowns: 
say Xy. 1. vanishes, we go over to a new basis by eliminating the 
unknown x, " from the old basis and introducing instead X, gy 
Going over to the new basis completes one step of the simplex 
method. 
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So far we have presented only some of the considerations 
poy for the understanding of the method of potentials ns 
‘4 the actual solution of the problem, it can be best explained 
y considering a specific example. 
Nun. We are given three source points, A;, 45, 45, with the 


a; = 60, a;—80, a= 100 
and four destination points, Bı, B», Bs, B4, having the requirements 
b; =40, bz = 60, b; =80, by =60 


The values cy (costs of transportation from A; to B, per ton of 
goods) are given by Table 6. 


Table 6 
1 2 3 4 
4 3 2 0 
0 2 2 1 


the finding of the first basic solution. 
corner rule leads to the result shown 
of the previous subsection). We denote 


i Solution. We begin with 
F this case the northwest 
Table 7 (cf. the example 
Table 7 


BSE: 


e function $ is as 


this solution by X for short. The value of th 
Ollows for it; — 


Lex, = 1 MERC ULLA EL 


42x04-2x 40+ 1 x 60— 420 


* £9 x40 r OORO 


To find the potentials it is necessary to solve the system of 
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equations 
a + By = ey1 = 1 
9, +B. =cy2=2 
95 + Ba =C22=3 
a2 + B3 =c23=2 
as + Bs =c33=2 
a3 + By =c34=1 


any solution of this system suiting us. As already noted, the value of 


one of the unknowns can be specified arbitrarily. Setting o; = 1, we 
find successively: 


Bi=0, Bo=1; a =2, B,=0, o,=2, B,——1 
We then compute the indirect costs Cig 
€73 =% + B3 = 1’ 
Ci 704 + By =0' 
c=% + By = 2’ 
€34,70; + By = 1' 
c31 =03 + B, = 2’ 
C32 = %3 + B2 = 3' 


It is more convenient- to carry out the above computation 
using Table 8 in which we first enter only the values of c, for 


ali the basic cells. Assuming x, = 1 and proceeding according to the 
above rule we compute the potentials v, and B, and, on entering 
them into the corresponding cells, we get Table 9. We then find the 


indirect costs c; entering them in the corresponding cells we get 
Table 10. 


Table 8 


[m 
N 
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Table 9 


Tu D 1 2A 
a 

1 1 2 

2 3 

2 1 


Table 10 


In fact, the three elementary stages we have considered could be 


P" oe ee the differences Sg = Spa Sot In this case 
gg=a-l= 2 
E = d 
S21 ELEC 2, 
Sag Oli = ails 
sna leg i a 
salare : 


5 ic unknowns 
Hence the function S expressed in terms oa 
assumes the form (8) 


5 = 420 + Dea + xia Di — X24 e 


are negative ones, the coefficient X31, 
Possible to decrease th 
Preserving the zero values © 
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We set x3, =p. Since the sums of the values of the unknowns in 
the rows and columns must remain unchanged, it is necessary to 
add p to or subtract it from the values of the basic unknowns 
in the following manner: 


Table 11 ` 


We had to compensate for the addition of p to x3, by subtracting 
P from x,, and in turn to compensate for this by the addition of. 
P to X,5 and so on until we came back to X31. , 

We stop here to make a remark concerning Table 11. Going in the 
table from cell to cell in that. sequence in which we compensate 
for p we get a closed broken line consisting of alternating horizontal 
and vertical links: the broken line is dotted in Table 11. One 
of the vertices of the broken line is in a nonbasic cell (x3, in 
this case) and the rest are in basic cells (not necessarily in all of them; 
thus there is no vertex in the cell x34 of Table 11). This broken 
line is called a Cycle or, to be more exact, a recalculation cycle 
corresponding to a given nonbasic cell. 

As is seen from Table 11, the condition that the unknowns should 
be nonnegative allows p to be increased only to 40. We set p= 40: 
three basic unknowns, Xii, X22 and x33, then vanish together. We 
choose one of them, X22, for example. So we introduce the unknown 
X31 into the basis and make the unknown x5, nonbasic. The new 
basic solution will be as follows: 


0 60 


40 0 60 


the value of the function S for it being 
420-2 x 40 = 340 
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(in (8) we set x3, = p — 40 and assume the values of the remaining 
unknowns to be zero). Going over to the new basis completes one 
Step in the procedure. We have succeeded in decreasing the value of 
the function S by 80 units. 

To sum up. by the beginning of the first step we had been 
Blven an initial basic solution X. The step itself consisted in going 


Over to new basic solution X and included the following stages: 


(1) finding the potentials a, B; and the indirect costs 
Em 
Gq = 0, +B, 
(2) computing the differences Spy = €, — Guy: 3 4 
(3) choosing a vacant cell corresponding to a negative difference 
Shq and constructing the recalculation cycle for this cell; 
_ (4) finding the new basic solution and the new value of the 
function s. y 
Notice that had all the differences Spa. proved nonnegative, that 
Would have meant that the basic solution x was optimal Then 


Stages (3), (4) would have become unnecessary. Thus the nonnegativity 


f all differences s, is the criterion for terminating the podm 
Let us return to the example under consideration, however. The 


next step begins with finding the potentials, indirect costs and comput- 
ing the differences s... The corresponding computation will look as 


Ollows: eg; 


0 

1 

v 

0 

t etta =n: 
543=3-3'=0 $22 =3 D 
ETE Lm Asa 
Sia 7 z PRONTA, 
So) = 4 ee 32 

+ X32 


S = 340 + 2X14 + 4x51 + 2X22 — X24 
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There are negative differences among Spg: S24 = — 1. We fix on 
the cell x;, and construct a recalculation cycle for it: 


40 O+p |60—p 


We conclude from the examination of the cycle that pcan be increased 
only to 60. Setting p = 60, we arrive at a new basic solution: 


3 20 60 


40 60 


The value of the function S is 


340 — 1 x 60 = 280 
for this solution. Going over from X to x completes the second 
step in the procedure, : 


The third step again begins with computing the potentials, indi- 
rect costs and differences Spy We have 


t3 
u 


Q 
N 
© 


o 
N 
Q 
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Now the differences s 


are all nonnegative. This means that the 


ra 


last basi i 
ic sol is i 
ution x is optimal and the corresponding (smallest) 


aa of the function S is 280. 
the optimal solution of the problem has been found: 


Ee = Sey = 
l 127 X347 X34 = 60, x33 — 20, x31 = 40; 


the remaining x; are zero. 


In c 1 ` 
onclusion we show once again, but now only schematically, 


t e 
he plan of solving the transportation problem: 


Finding the first basic solution 


| 


Finding the potentials o, Pı 


m 


Computing the indirect costs 6 ng 


Computing the differences 


ert 
Sp = Sra ~ pa 
(the coefficients of nonbasic unknowns in 


the expression for Sy 


Choosing a negative difference s, and 
constructing the recalculation cycle for 
the corresponding vacant cell 


I 


Finding a new basis 


is u Imi oe m 


E If the differences Sp are all nonnegative, t 
Cn is optimal. 
. Justifying the method of potentials. The formula 


=f iei d 
Spa = Cpa Cpa 


awe the function s 
<nowns is the pivot o 
Meee this formula requires a deeper stu 

le and we have agreed to discuss it se 


hen the last basic 


(7) 


to be expressed in 0 
d described above. 
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it now. To that end we return to the study of Table 11 
Table 11 


which contains the recalculation cycle for the cell X31. We used 
this table to find p and thereby go over to a new basic solution. 
But the table allows still another problem to be solved: it helps 
to determine what are the coefficients of the unknown X3, in the 
expressions for the basic unknowns. 

We first of all make the following remark. If the values of nonbasic 
unknowns are indicated, then this unambiguously determines the 
values of the basic unknowns. But it is just Table 11 that gives 
an idea of what happens to the values of basic unknowns when 
the nonbasic unknown X3; is assigned a value p (while assuming 
the remaining nonbasic unknowns to be zero). The study of the table 


shows that *31 enters the expressions for the basic unknowns having 
the following coefficients: 


Xni 
aioe Met EI 
cubus (9) 
ar lotos. 


=A S 


In the remaining basic unknowns (namely, in X34) Xa, has the 
coefficient 0. We see that when making the round of the cycle the 
coefficients involved alternatively assume the values th=]. 

It now remains to answer the question: What coefficient has 
X3, in the expression for the function S in terms of nonbasic 
unknowns? Well, this expression 


S = 420 + 513x134 S14X14 + 554x514 524X34 + 
T $3131 + 535x3; 


is produced when the basic unknowns in the original expression 
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for S 

S=Zex, (i=1, 2, 3; j=1, 2, 3, 4) 
E substituted for by their expressions in terms of nonbasic ones. 
n the given case it is necessary that x11, X12; X22, X23, X33, X34 
Should be substituted for by their expressions in terms of x13, X14 
ee X24, X31, X32. In virtue of (9) the total coefficient of x3, will 


$31 = C31 — Cii + C12 — C22 + C23 — C33 


Hence 
S31 = c31 — (04 + Bi) + (1 + B2) — (2 + Ba) + 


+ (at; + B3) — (os + B3) = C31 — (@3 + Bi) 
or 
Si 63 Lae [T 


- But that is just formula (7) (for the given case). The reasoning that 
has led us to the last equality is sufficiently demonstrative, and 
there is no need to repeat it for the general case. ) 
We conclude by saying a few words about other methods of solving 
the transportation problem. The method of potentials discussed 
above tends to be used in hand computation since one of its stages, 
Namely that of constructing the recalculation cycle, is difficult to 
Tealize in the computer. With a small number of suppliers and 
Consumers (within some ten) this method leads to an optimal solution 
in a reasonable time. If, however, the number of suppliers and that 
9f consumers are large enough (and this is what happens in appli- 
Cations, where these numbers are of an order of hundreds), resort- 
ing to a computer is inevitable. Algorithms computer programmes 


are based on considerations other than those underlying the method 
der can find more about these 


of potenti i d rea 

als. The interested re [ 
algorithms in special literature devoted to the transportation 
Problem, 
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METHOD OF SUCCESSIVE APPROXIMATIONS 


N. VILENKIN, D.SC. 


This book explains in a popular form the methods of appro- 
ximation. solutions of algebraic, trigonometric, model and 
other equations. 

Intended for senior schoolchildren, polytechnic students, ma- 
thematics teachers and for those who encounter solutions of 
equations in their practical work. In the course of exposition 
some elementary ideas about higher mathematics are 
introduced in the book. About 20 solved exercises are 
included in the appendix. 

The book has already been translated into Spanish and 
French. 


BASIC COMPUTATIONAL MATHEMATICS 


V. DYACHENKO, D.SC. 


The present text covers the most eleme 
ideas underlying modern computational 
problems in mechanics and mathematical 
with the construction and investigation of 
tational algorithms, 
The presentaiion is simple and does not require that the 
reader have a substantial mathematical background. The text 
is designed for students of natural-science departments of colle- 
ges and universities. It will also be of interest to a broad range 
of physicists and engineers serving as an introductory course 
in computational mathematics. 


ntary concepts and 
methods for solving 
physics. It also deals 
appropriate compu- 
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FUNDAMENTAL THEOREM OF ARITHMETIC 


L. KALUZHNIN, D.SC. 


This booklet is devoted to one of the fundamental proposi- 
tions of the arithmeticof rational whole numbers —the 
unique factorisation of whole numbers into prime multipliers. 
It gives a rigorous and complete proof of this basic fact. It 
is shown that uniqueness of factorisation also exists in arithme- 
tic of complex (Gaussian) whole numbers. The link between 
arithmetic of Gaussian numbers and the problem of represent- 
ing whole numbers as sum of squares js indicated. An example 
of arithmetic in which uniqueness of expansion into prime 
multipliers does not hold is given. The booklet is intended for 
senior schoolchildren. It will help to acquaint them with the 
elements of number theory. It may also be useful for secondary 
school teachers. 


INDUCTION IN GEOMETRY 


L.GOLOVINA, D.SC. and I. YAGLOM, D.SC. 


This booklet deals with various applications of the method 
of mathematical induction to solving geometric problems and 
was planned by the authors as a natural continuation of I. S. So- 


minsky's booklet “The Method of Mathematical Induction" ` 


published by Mir Publishers in 1975. It contains 37 worked 
examples and 40 problems accompanied by brief hints. 


TER OU ONCE 
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The book tells about the relation of systems of linear inequali- 
ties to convex polyhedra, gives a description of the set of all 
solutions of a system of linear inequalities, analyses the ques- 
tions of compatibility and incompatibility; finally, it gives an 
insight into linear programming as one of the topics in the 
theory of systems of linear inequalities. The last section but one 
gives a proof of the duality theorem of linear programming. 
The book is intended for senior pupils and all amateur 
mathematicians. 


Mir Publishers 
Moscow 


